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ABSTRACT 


This  report  is  a  review  of  the  literature  concerning  the  excitation 
of  structures  by  random  pressure  fluctuations  in  a  turbulent  bound¬ 
ary  layer  and  the  consequent  acoustic  radiation .  Such  radiation  is 
termed  flow  noise.  The  purpose  of  the  report  is  to  determine  those 
aspects  of  analysis  and  experiment  on  the  subject  that  are  pertinent 
to  the  generation  of  self  noise  in  a  sonar .  The  report  is  concerned 
with  the  physical  mechanisms  governing  flow  noise,  with  the  pre¬ 
diction  of  the  magnitude  and  distribution  in  frequency  of  flow  noise, 
and  with  the  discovery  of  measures  that  may  be  taken  to  reduce  its 
influence  upon  sonar  self  noise .  The  outstanding  research  problems 
remaining  in  the  field  are  pointed  out,  and  recommendations  are 
made  for  the  direction  of  further  research . 


vii 


3riliur  2Ll.ittlr.ilnr. 


I.  INTRODUCTION 


Sonar  self  noise  is  a  combination  of  many  effects  in  addition  to  flow 
noise.  These  include  internal  electrical  noise  in  the  sonar  system,  noise 
from  own  ship's  machinery,  sea  ambient  noise,  discrete  tones  from  singing 
propellers  and  appendages,  general  rattles  and  bangs  from  loose  equipment, 
and,  in  the  case  of  surface  ships,  cavitation  noise  from  the  propellers  or 
from  cavitation  or  air  bubble  impingement  on  the  sonar  dome.  We  do  not 
concern  ourselves  with  any  of  these  additional  effects  in  this  report .  Our 
interest  is  directed  specifically  to  noise  resulting  from  pressure  fluctuations 
in  the  turbulent  boundary  layer  around  the  sonar  dome  at  high  ship  speeds . 

Such  noise  contributes  significantly  to  sonar  self  noise  on  deeply  submerged 
submarines  traveling  at  high  speed  and  plays  an  important  role  in  surface  ship, 
torpedo,  and  towed  sonar  self  noise  as  well. 

The  literature  relevant  to  our  problem  covers  a  period  of  about  fifteen 
years.  Most  of  it  is  directed  toward  corresponding  problems  in  aerodynamic 
sound  involving  skin  vibration  and  internal  vibration  and  sound  in  missiles  and 
aircraft  fuselages.  Although  much  of  the  work  on  aerodynamic  sound  is  di¬ 
rectly  applicable  to  the  sonar  self  noise  problem  rhere  are  significant  differ¬ 
ences  in  the  sonar  problem  because  of  the  relatively  low  flow  speeds  involved 
and  the  high  density  of  the  water  medium .  One  of  our  purposes  is  to  point  out 
the  effect  of  these  differences . 

Sections  II,  III,  and  IV  of  this  report  briefly  survey  background  ma¬ 
terial  pertinent  to  the  principal  subject.  Section  II  is  concerned  with  semi- 
empirical  studies  of  mean  velocity  distributions  in  a  turbulent  boundary  layer; 
these  studies  were  originally  performed  to  determine  frictional  drag.  In 
Section  III  a  framework  for  the  statistical  analysis  of  turbulent  boundary  layer 
flows  is  constructed,  and  in  Section  IV  the  foundations  of  the  theory  of  aero¬ 
dynamic  sound  are  discussed . 

From  these  foundations  certain  analytical  relations  between  wall 
pressure  fluctuations  and  velocity  fluctuations  can  be  established  and  are  dis¬ 
cussed  in  Section  V.  Because  of  formidable  difficulties  that  stem  from  the 
non-linearity  of  the  relations,  only  crude  results  can  be  obtained.  Of  much 
greater  significance  are  the  actual  measurements  of  wall  pressure  fluctuations 
reviewed  in  Section  VI.  These  provide  the  statistics  of  the  wall  pressure  field 
that  is  the  random  forcing  function  of  our  problem .  The  results  in  this  section 
are  applicable  to  flush-mounted  sonar  transducers  where  the  wall  pressure 
field,  termed  pseudo-sound,  contributes  directly  to  transducer  self  noise. 
Unfortunately ,  the  engineering  usefulness  of  experimental  results  obtained  with 
transducers  of  finite  size  has  been  clouded  by  controversies  over  their  inter¬ 
pretation  . 
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Analytical  studies  of  the  response  of  infinite  plates  and  periodically 
simply  supported  plates  to  wall  pressure  fluctuations  and  of  the  consequent 
acoustic  radiation  from  the  plates  are  treated  in  Section  VII .  Section  VIII 
discusses  the  results  of  measurements  taken  to  test  the  validity  of  these  an¬ 
alyses.  All  tests  reported  were  conducted  by  installing  thin  plates  in  the  walls 
of  wind  tunnels  and  measuring  plate  vibrations  and  acoustic  radiation  fields  as 
functions  of  tunnel  velocity  and  plate  dimensions .  Several  of  the  analyses  and 
experiments  included  studies  of  the  effect  of  structural  damping  upon  plate 
vibration  and  acoustic  radiation.  The  results  of  both  the  analyses  and  the 
experiments  are  summarized  in  Section  IX.  Section  X  gives  a  summary  and 
conclusions  related  to  the  report  as  a  whole  and  includes  recommendations 
for  further  research  on  aspects  of  flow  noise  related  to  sonar  self  noise , 

For  clarity  of  presentation,  notations  used  by  the  original  investi¬ 
gators  have  been  made  to  conform  to  a  single  notation  scheme .  This  is 
particularly  true  in  the  identification  of  pressure  correlation  and  spectral 
functions ,  as  no  uniformity  has  evolved  in  the  literature .  In  some  of  the 
analyses ,  the  original  derivations  have  also  been  altered  in  an  effort  to  lend 
coherence  to  the  report,  but  care  has  been  taken  to  retain  the  basic  assump¬ 
tions  and  conclusions  of  the  originators.  In  particular.  Section  VII,  the 
analysis  of  Ffowcs  Williams  and  Lyon  (.1963)  has  been  recast  in  the  framework 
originally  derived  by  Kraichnan  (1957);  this  is  done  in  order  to  explain  the 
marked  difference  in  the  dependence  of  radiated  acoustic  power  per  unit  area 
upon  displacement  boundary  layer  thickness  predicted  in  these  papers.  How¬ 
ever,  the  interpretations  and  comparisons  of  results,  together  with  the  final 
conclusions  and  recommendations,  are  the  responsibility  of  the  present  author. 

The  author  wishes  to  acknowledge  the  help  of  a  number  of  colleagues 
whose  comments  aided  the  preparation  of  this  report.  Remarks  by  J.  E. 
Barger,  K.  L.  Chandiramani,  G.  M.  Corcos,  I.  Dyer,  J.  E.  Ffowcs  Williams, 
F.  J.  Jackson,  R.  H.  Lyon,  G.  Maidanik  and  M.  Strasberg  were  of  particular 
value . 
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II.  SEMI  -EMPIRICAL  RELATIONS  FOR  MEAN  VELOCITY 
DISTRIBUTION  IN  A  BOUNDARY  LAYER 


The  foundations  of  flow  noise  research,  although  not  recognized  as 
such  at  the  time,  were  laid  in  the  early  investigations  of  turbulent  boundary 
layer  flow  by  Prandtl,  von  Kdrman,  Nikuradse,  and  Reichart.  The  work  was 
semi -empirical  in  nature  and  directed  toward  determining  the  frictional  resis¬ 
tance  to  turbulent  flow  in  pipes  and  over  flat  plates .  Later  investigations 
treated  the  effects  of  surface  roughness  and  flows  over  various  body  shapes 
where  significant  pressure  gradients  occur.  A  detailed  presentation  of  the 
results  is  given  in  the  book  by  Schlichting  (1955).  We  summarize  the  results 
of  primary  interest  here . 

It  is  now  well  established  (Coles  1956)  that  an  intermediate  region 
of  the  turbulent  boundary  layer  produced  by  flow  over  both  smooth  and  rough 
surfaces  satisfies  a  uniform  velocity  distribution  relation  known  as  the  "law 
of  the  wall" : 


u(y)  =  I 

v*  * 


In  (  y  — )  +  C 

v 


(1) 


where  u(y)  is  the  time -averaged  tangential  velocity,  y  is  the  distance  normal 
to  the  wall,  and  v  is  the  kinematic  viscosity  of  the  fluid.  The  constants 
K  =^0.40  and  C  *=5.1  are  found  to  be  valid  for  all  boundary  layer  flows.  The 
quantity  v*  is  called  the  friction  velocity  and  is  defined  as 


} 


1 

I1 

fi 

11 

I 


v 


(2) 


where  tq  is  the  shear  stress  at  the  wall  and  pQ  is  the  fluid  density,  a  constant 

v, 

for  low  Mach  number  flows .  The  quantity  y  —  is  evidently  a  Reynolds  number; 

Eq.  (I)  is  valid  when  it  exceeds  70.  If  we  define  the  boundary  layer  thickness 

6  as  that  value  of  y  for  which  u  =  0.99U,  where  U  is  the  constant  free -stream 
velocity,  Eq.  (1)  can  be  written  as 

T.  =  5-75  ‘°8  10  (6>  +  ^  (3) 
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When  y  ~  <  5,  a  laminar  sublayer  is  found  for  which  ~  =  V  ~  . 
We  denote  its  thickness  6  by 


.  ,  v 

f  v* 


If  the  surface  rouglmess  is  less  than  6^.  and  the  wall  is  essentially  smooth,  we 
may  approximate  Eq.  (3)  by  the  relation 


—  =  8.74  (y  WP 


under  the  assumption  that 


Two  other  thicknesses  are  important;  one  is  the  displacement  thick¬ 
ness  of  6*  previously  referred  to: 


j  (U-u)  dy , 


and  the  other  is  the  momentum  thickness  9: 


e  =  U5  j  u(u'u>  dv- 

o 


6  7 

Using  Eq.  (6)  it  follows  directly  that  6*  =  — •  and  fi  =  —  6. 

o  12 


In  the  absence  of  a  pressure  gradient  along  the  wall. 


2  dP 


t  (x)  =  pU  -7- 
o  dx 


where  x  is  the  coordinate  along  the  wall  in  the  direction  of  flow. 


Equation  (5)  defines  v*  empirically  for  a  given  u(y),  and  we  can 
obtain  successively  the  results 


(10) 


.0225  (j 


v_\l/4 

U&/ 


6<x)  = 


0.37  x 


(11) 


and 


6<x) 


.036  x 


()2) 


The  dimensionless  friction  coefficient  C 


f 


I 

o 


1  p  u2 
—  o 


becomes 


Cf  =  .0592  Q-y1'3  (13) 

with  a  minor  experimental  adjustment  of  the  coefficient.  The  Reynolds  number 
Ux 

—  increases  with  distance  x  from  the  leading  edge  of  the  wall .  Here  we  have 

assumed  the  boundary  layer  tripped  at  the  leading  edge  so  that  it  is  Turbulent 
along  its  length . 

In  point  of  fact,  the  1/7-power  velocity  distribution  law,  Eq.  (5),  is 
valid  only  for  smooth  plates  at  moderate  Reynolds  numbers .  We  have  used  it 
for  simplicity  to  permit  rapid, approximate  calculations  of  turbulent  boundary 
layer  thicknesses  and  resistance  coefficients .  More  accurate  calculations  are 
made  in  the  references  cited  from  the  more  general  relation  of  Eq .  (3) .  The 
friction  velocity  v*  is  the  critical  parameter.  It  is  dependent  upon  wall  rough¬ 
ness  and  must  be  obtained  experimentally,  either  from  direct  measurement  of 
wall  shear  stress  tq  or,  as  is  more  common,  from  momentum  balance  calcu¬ 
lations  based  upon  measured  values  of  mean  velocity  profiles  and  pressure 
gradients  along  the  wall . 

It  is  physically  evident  that  the  wall  shear  stress  and  friction  coef¬ 
ficient  are  markedly  increased  if  the  roughness  exceeds  the  thickness  of  the 
laminar  sublayer  for  a  smooth  wall.  Using  Eqs.  (2),  (4)  and  (13)  rhe  approx¬ 
imate  relation 


^  29.2 


^Ux  V/10 

V7/ 


<H) 


v 


is  obtained  as  a  criterion  for  admissible  roughness.  The  right-hand  side  is  a 

Ux  .5 

slowly  varying  function  of  Reynolds  number.  For  R  =  —  between  10  and 

10  it  is  roughly  equal  to  10C.  Note,  however,  that  this  is  a  resistance  or 
drag  criterion.  As  yet  it  says  nothing  directly  about  the  effect  of  roughness 
on  turbulent  boundary  layer  pressure  fluctuations . 

Coles  (1956)  proposed  an  extension  of  the  uniform  velocity  distribu¬ 
tion  law  of  Eq.  (1)  from  an  extensive  study  of  turbulent  boundary  layer  measure¬ 
ments  involving  mean  velocity  and  pressure  gradients ,  smooth  and  rough  sur  - 
faces,  and  separated  and  re -attached  boundary  layers.  His  semi -empirical 
expression,  known  as  the  'law  of  the  wake,"  is 


u 

v* 


l 

K 


In 


(??> 


(15) 


The  function  w(y/6)  is  the  additive  wake  effect  for  the  outer  portion 
of  the  boundary  layer.  Coles  tabulates  values  for  w  obtained  from  fitting  exper¬ 
imental  data;  however,  a  close  approximation  is 


w  (|)  «  1  +  sinn(y/6  -  1/2). 


<16) 


The  constants  K  and  C  are  the  same  as  in  Eq.  (1).  The  parameter  n 
is  in  general  a  function  of  the  streamwise  distance  but  is  constant  (^  0.55) 
for  a  constant  mean  pressure  flow .  It  is  given  by  a  simple  transcendental 
equation 

211  -  ln(l  -HI)  =  K  —  -  In  Q—  J  -  KC  -  In  K.  ^ 

The  boundary  layer  thickness  is  redefined  implicitly  in  Cole's  form¬ 
ulation  by  the  requirement  that 

I  GO dw  =  1 

o 

used  to  normalize  the  wake  function  w .  It  follows  from  the  definition  of  6* 
of  Eq.  (7)  and  Eq.  (15)  that 


K 


6*U 

6v* 


i  +  n 


(18) 
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Eq.  (15)  is  a  valid  representation  of  the  velocity  profile  everywhere 

AV  1  f/\  <■!«/>  1  r\  •  •«  rt  m  /*  •  aUl  Ilf  '  *-l»  *  —  4-t.  —  Z  -»«.  L>  1  O  ...t.  «.! 1  *  —  —  _ .  __ 

vAiviuui  cvj  tut  lcinixxiax  cuuiayci  ,  vt  ilim.ii  uic  latuiiiai  duuxayci  }  uur  linear 

relation 


u 

v* 


is  assumed  to  hold.  This  relation  is  merely  an  expression  of  shear  stress 
proportional  to  a  linear  velocity  gradient  valid  for  simple  laminar  viscous  flow. 

Cole’s  law  of  the  wake  is  currently  in  successful  use  for  determining 
mean  properties  incidental  to  boundary  layer  noise  studies .  All  five  of  the 
dimensionless  parameters 


_U 

v* 


6  *  fiv*  6*  U 
6  ’  v  ’  v 


and  II 


can  be  determined  if  any  two  are  known.  For  example,  suppose  that  a  velocity 
profile  u(y)  is  measured  in  a  particular  boundary  layer  experiment.  Then,  by 
curve  fitting  the  logarithmic  part  of  Eq.  (15),  the  friction  velocity  v*  can  be 
determined.  The  displacement  boundary  layer  thickness  can  be  computed 
from  u(y)  and  its  maximum  value  U.  With  y  ~  6  and  u  =  U  in  Eq.  (15)  the 
parameter II  can  be  determined.  Then  from  Eq.  (18)  the  boundary  layer  thick¬ 
ness  6  can  be  computed . 

In  pipe  or  duct  flow ,  v*  is  frequently  measured  directly  in  terms  of 
the  static  pressure  drop  over  a  length  of  uniform  test  section .  Assuming  that 
this  gradient  is  small,  the  use  of  n^O.SS  and  6  equal  to  the  pipe  radius  per¬ 
mits  one  to  treat  Eq.  (15)  as  a  close  approximation  to  the  actual  velocity  dis¬ 
tribution.  But  more  important,  perhaps,  is  the  fact  that  Cole’s  law  of  the  wake 
is  redundant;  i.e.,  more  measurements  are  practical  than  are  actually  required. 
This  has  permitted  the  establishment  of  the  important  experimental  fact  that 
fully  developed  turbulent  pipe  or  duct  flows  satisfy  the  same  mean  velocity  dis¬ 
tribution  law  as  external  turbulent  boundary  layer  flows .  By  inference  it  is 
highly  plausible  that  their  fluctuating  velocity  components  behave  similarly. 
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III.  THE  STATISTICAL  ANALYSIS  OF  TURBULENT  BOUNDARY  LAYER  FLOWS 


C  lly  one  statistical  property  of  a  turbulent  flow  has  been  introduced 
thus  far,  namely,  the  time  average  of  velocity.  G.I.  Taylor  (1935)  (1936)  was 
the  first  to  apply  the  methods  of  statistical  mechanics  to  the  study  of  turbulence. 
Higher  order  ensemble  averages  of  velocity  components  were  studied  within  the 
framework  imposed  by  the  laws  of  conservation  of  mass  and  momentum  for  in¬ 
compressible  flow.  The  theory  rests  upon  two  basic  hypotheses:  first,  the 
ergodic  hypothesis  that  an  ensemble  average  is  equivalent  to  a  spatial  average; 
and,  second,  Taylor's  hypothesis  that  a  time  average  at  a  fixed  point  in  a  turbu¬ 
lent  flow  is  equivalent  to  a  spatial  average  taken  along  a  line  in  the  direction  of 
mean  flow.  Many  measurements  of  velocity  fluctuation  with  time  have  been  made 
using  hot  wire  techniques  as  necessary  supplements  to  the  theory. 

The  books  by  Batchelor  (1953),  Townsend  (1956),  and  Hinze  (1959) 
give  detailed  presentations  of  both  theoretical  and  experimental  results .  Although 
they  contain  little  theory  and  no  experiments  at  all  related  to  pressure  fluctua¬ 
tions  in  a  boundary  layer,  they  do  set  the  foundations  on  which  our  understanding 
of  flow  noise  has  developed.  A  brief  statement  of  features  of  this  work  pertinent 
to  flow  noise  research  is  therefore  necessary. 

A  velocity  correlation  tensor  is  defined  as 


R..(r )  =  u.(x  )u.  (x  +  r ) 

ij  i  j 


(19) 


where  uj(x  )  is  a  velocity  component  at  the  point  x.  The  bar  denotes  ensemble 
average,  e.g.. 


Ria(r  )  =  f  ux(x)u2(x  +  r )  P(u1,u2)  dUj  du2 

where  P(ux  ,u2)  is  a  joint  probability  density  function.  Equation  (19)  states  that 
Rjj  is  independent  of  x,  implying  that  the  flow  is  spatially  homogeneous .  T.ie 
ergodic  hypothesis  implies  that 


R..(r)  = 


lim 


I'- 


i(x)u.(x+F) 


dx 


where  V  is  the  volume  of  integration  in  x  space . 
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The  Fourier  transform 


Mkl  = 

1J  '  m 


1  In  /~*\  -  ik ■  r 

—  /  R..(r)  e 

nf  J  « 


dr 


(20) 


is  called  an  energy  spectrum  tensor.  It  has  the  form  of  a  wavenumber  spectral 
density,  for  by  taking  the  inverse  transform  we  obtain 


V°>  * 


■/v‘ 


k )  dk 


(21) 


Alternatively,  we  can  take  the  Fourier  transform 


u  (k)  =  —  ju 
(2tt)3  J 


—  - ik  •  X  — 

(x) e  dx 


(22)f 


where  we  assume  u.(x)=0  exterior  to  a  volume  V.  It  follows  that 

j 


Mk)  =  lim  u*  (k)u.(k) 

y  \r_»— .  V  l  J 


(23)' 


where  u  *  is  the  complex  conjugate  of  uj.  Those  statements  can  be  generalized 
directly  to  include  time  variable  t  and  its  corresponding  frequency  variable  tr, 
or  specialized  to  reduce  to  two  space  variables .  frequently  the  quantities  are 
normalized,  e.g.,  by  dividing  Rlg  (r  )  by  (u3)*  (u*)^ . 

It  is  helpful  to  keep  in  mind  that  the  physical  properties  of  the  fluid 
are  the  random  variables;  hence,  spatial  and  time  derivatives  of  velocity  yield 
space  and  time  derivatives  of  the  velocity  correlation,  but  only  k  and  w  multiples 
of  the  spectrum  tensor .  For  this  reason,  spectrum  relations  evolving  from  sta  - 
tistical  treatment  of  momentum  and  continuity  equations  are  sometimes  more 
readily  treated  analytically. 


t  Equations  (22)  and  (23)  are  formal  statements,  since  the  integral  in  Eq.  (22) 
diverges  in  general  for  real  velocity  distributions .  A  more  rigorous  approach 
is  given  by  Batchelor  (1953). 


9 


3rtI|ur  Sl^UttlcJhtc. 


As  an  example,  we  apply  the  above  concepts  to  the  fluctuating  pres¬ 
sure  p(x,t).  If  the  pressure  is  statistically  homogeneous  in  space  and  time,  the 
pressure  correlation  is 


p(x,t)  p(x+r,t+T)  =  R(r,T) 


(24) 


and  R(0, 0)  is  the  mean  square  pressure.  The  Fourier  transform 


Hk)  = 


1 


<2n)3 


flKr-.Oe^ 


dr 


(25) 


is  called  the  wavenumber  spectral  density  of  the  pressure;  the  Fourier  transform 

m  =  4r  m°.  t>  e '  il"Tdf  <26> 


is  called  the  frequency  spectral  density  of  the  pressure;  and  the  Fourier  transform 

H  £>)  =  Jr(  t ,  t)  e"lUJT  dT  (27) 

is  called  the  frequency  cross  spectral  density  of  the  pressure.  In  dealing  with 
the  fluctuating  pressure  exerted  by  a  turbulent  boundary  layer  on  a  flat  wall  (wall 
pressure)  we  shall  customarily  restrict  vectors  x  and  r  to  lie  in  the  plane  of 
the  wall  by  setting  their  normal  components  x2  and  r2  equal  to  zero .  The  factor 
(2n)-3  in  Eq.  (25)  is  then  replaced  by  (2n)-2  and  the  integration  is  carried  out 
over  rx,  r3  only.  The  vector  k  has  the  components  ka ,  k^  only  in  this  case. 

We  must  bear  in  mind  that  it  is  practical  to  measure  only  time  varia¬ 
tions  of  physical  quantities  such  as  velocity  and  pressure.  Yet  all  of  the  aver¬ 
ages  discussed  above  are  ensemble  averages .  Clearly  the  physical  circumstances 
of  a  turbulent  flow  must  be  such  that  both  an  ergodic  hypothesis  and  Taylor’s  hy¬ 
pothesis  can  be  invoked  if  there  is  to  be  any  relationship  between  experiment  and 
the  statistical  theory  of  turbulence.  This  point  is  fundamental  to  boundary  layer 
noise  research  as  well  as  to  the  general  study  of  turbulence.  We  need  to  examine 
the  applicability  of  these  hypotheses  to  fluctuating  velocities  in  order  to  gain  in¬ 
sight  into  their  applicability  to  fluctuating  pressures . 
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Let  us  consider  first  the  generation  of  turbulence  by  placing  a  fine 
mesh  grid  across  the  upstream  side  of  the  test  section  of  a  wind  tunnel.  Let  U 
be  the  mean  tunnel  flow  velocity  and  consider  it  directed  in  the  positive  xx  direc¬ 
tion.  Let  uj(x,t)  be  the  fluctuating  turbulent  velocity  components .  If  the  tunnel 
has  been  in  steady  operation  for  a  reasonable  length  of  time,  it  is  clear  that  a 
velocity  correlation  such  as 


R..(x,x  +  r,T)  =  u.(x,t)  u.(x  +  r*,t  + 1) 

is  not  dependent  upon  t.  Hence  we  may  make  the  ergodic  hypothesis  that  Rjj  is 
equal  to  the  time  average 


T 


u.(x,t)  u.(x  +  r,t  +  t)  dt 


taken  over  a  sufficiently  long  time  T.  This  time  average  can  be  measured  by 
placing  hot  wires  at  the  points  x  and  x  +  r,  and  performing  the  time  delay  and 
integration  operations  on  the  output  signals . 

Let  us  now  consider  a  slowly  decaying  (or  varying)  turbulent  velocity 
field  being  converted  at  velocity  Uc  in  the  positive  Xj direction.  We  have,  roughly, 
that 


u.(xi,x2,x3,t) 


Strict  equality  would  imply  that  the  turbulent  velocities  are  frozen  in  time  in  a 
reference  frame  moving  with  the  convection  velocity  Uc.  Taylor’s  hypothesis 
is  that  this  relation  is  valid  when  applied  to  velocity  correlations  and  that,  fur¬ 
ther,  the  turbulence  is  spatially  homogeneous  in  the  moving  frame  of  reference. 
Thus  for  x2  =  x3  =  r2  =  r5  =  0,  say. 


u.(xl,t)uj(x1 


+  rx,t  +  T) 


-ucT,°) 


(28) 
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We  now  have  a  means  of  determining  experimentally  the  correlations 
required  in  the  statistical  analyses.  Favre,  Gaviglio,  and  Dumas  (1953)  demon¬ 
strated  the  validity  of  Eq.  (28)  for  grid  turbulence  by  showing  that  peak  values 
of  the  velocity  correlations  were  obtained  for  rx  =  Uct  .  This  form  of  turbulence 
is  essentially  isotropic  in  the  sense  that  the  velocity  correlation  Rij(r)  is  invariant 
under  rigid  rotations  and  reflections  of  ui(x),  Uj(x-i  r)  and?.  Physically,  this 
means  that  there  is  no  directional  preference  for  turbulent  velocities  in  free  tur¬ 
bulence  produced  by  grids .  Lin  (1953)  showed  analytically  that  Taylor’s  hypothesis 
is  valid  for  isotropic  turbulence,  provided  the  root-mean -square  velocity  fluctua¬ 
tions  are  much  less  than  the  convection  velocity  Uc . 

For  turbulent  shear  flows,  such  as  occur  in  the  boundary  layer,  there 
is  less  theoretical  justification  for  applying  Taylor's  hypothesis  to  velocity  cor¬ 
relations.  However,  Favre,  Gaviglio,  and  Dumas  (1958)  demonstrated  that  it 
does  apply  if  the  convection  velocity  Uc  is  taken  equal  to  the  mean  flow  velocity 
at  the  distance  x2  normal  to  the  wall  at  which  the  hot  wire  measurements  were 
taken.  In  boundary  layers,  no  assumption  of  general  spatial  homogeneity  can  be 
made.  However,  it  is  customary  to  assume  that  the  turbulent  boundary  layer  is 
spatially  homogeneous  in  planes  parallel  to  the  wall .  This  assumption  is  made 
physically  tenable  when  one  considers  that  similarity  must  persist  in  the  bound¬ 
ary  layer  over  lengths  equal  to  many  boundary  layer  thicknesses . 

We  must  emphasize  that  none  of  the  extensive  studies  of  velocity  cor¬ 
relations  in  turbulent  boundary  layers  give  us  any  immediately  useful  information 
about  the  pressure  fluctuations .  They  do  serve  two  purposes,  however.  First, 
they  set  a  pattern  of  analysis  and  experiment  on  which  to  model  studies  of  pres  - 
sure  fluctuations.  Second,  as  we  shall  see,  pressure  fluctuations  are  related 
to  velocity  fluctuations  through  the  dynamical  equations  of  motion;  hence,  informa¬ 
tion  that  may  be  inaccessible  through  pressure  measurements  can,  in  principle, 
be  sought  from  velocity  measurements. 
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IV.  THE  FOUNDATIONS  OF  THE  THEORY  OF  AERODYNAMIC  SOUND 


Lighthill  (1952)  published  a  classic  paper  giving  the  theory  of  acoustic 
radiation  from  free  turbulence.  This  paper  is  the  starting  point  for  all  subse¬ 
quent  investigations  of  flow  noise.  In  it  he  made  the  first  quantitative  estimates 
of  the  intensity  of  acoustic  radiation  from  free  turbulence,  showing  it  proportional 
to  the  eighth  power  of  the  mean  flow  Mach  number  for  subsonic  flows.  He  ex¬ 
pressed  the  equations  of  motion  in  the  form 


_^!p. 
3y  d\ 
i  '  i 


_L 

c2  St2 


2 

c 

o 


d2T 
_ y_ 

dx.  dx. 
1  J 


(29) 


where  P  is  the  fluid  density  and  cG  is  the  speed  of  sound  in  the  quiescent  portion 
of  the  fluid.  The  summation  convention  is  used  where  the  same  index  appears 
twice  in  a  term,  e.g. 


*2P 


dx.  ox. 

i  i 


a2p  d2p  a2p 


The  tensor  T. .  is  defined  as 
ij 


T..  =  Pu.u.  -  7..  -  C2P6..  (30) 

ij  i  j  ij  o  ij 


where  Uj  is  a  velocity  vector,  6jj  is  the  Kronecker  delta  (=  1  for  i  =  j  and  =  0 
for  i  f  j).  The  stress  in  the  fluid  is 


=  -  P 


y 


ij 


where  p  is  the  pressure  and  7j  represents  the  viscous  components  of  stress 
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In  a  region  of  turbulence  the  dominant  term  of  Ty  is  the  Reynolds 
stress  PujUj.  In  the  quiescent  portion  of  the  fluid  Ty  may  be  assumed  to  vanish. 
Thus  Eq.  (29)  can  be  interpreted  as  a  wave  equation  for  the  radiation  of  density 
fluctuations  outward  from  a  distribution  of  quadrupoles  Ty  representing  the 
turbulence.  By  the  application  of  Kirchhoff’s  formula  for  a  retarded  potential 
to  Eq.  (29),  Lighthiil  demonstrated  that  quadrupole  or  free -turbulence  radiation 
intensity  obeys  the  above  law  involving  the  eighth  power  of  the  Mach  number . 

On  one  hand,  his  result  is  of  considerable  importance  in  the  study  of 
aircraft  jet  engine  noise.  On  the  other  hand,  it  serves  to  show  negatively  that 
free  turbulence  noise  is  unimportant  in  underwater  applications,  where  the  sound 
velocity  is  high  and  vehicle  speeds  are  low. 

Curie  (1955)  extended  Lighthiil ’s  theory  to  include  the  effect  of  a  rigid 
body  adjacent  to  a  region  of  turbulence.  He  showed  that  the  intensity  of  radiation 
from  surface  pressure  fluctuations  on  the  body  induced  by  the  turbulence  should 
be  proportional  to  the  sixth  power  of  the  mean  flow  Mach  number.  The  dominant 
contributor  to  radiated  sound  in  this  case  is  a  dipole  distribution  at  the  body  sur¬ 
face. 


Curie's  theory  was  applied  by  Phillips  (1956a)  to  predict  the  intensity 
of  aeolian  tones  created  when  a  body  in  a  flow  field  sheds  a  regular  vortex  wake . 
The  radiation  from  cylinders  whirled  in  still  air  or  placed  in  subsonic  wind  tun¬ 
nels  was  measured  by  Phillips  (1956a),  Gerrard  (1955),  and  Keefe  (1961);  their 
findings  generally  confirm  the  law  concerning  the  sixth  power  of  the  Mach  number, 
provided  proper  account  is  taken  of  the  length  of  the  cylinder  over  which  the  shed 
vortices  are  correlated  in  phase.  Dipole  radiation  of  this  type  can  be  important 
in  underwater  applications .  Ross  (1964)  showed  that  it  is  probably  the  principal 
source  of  sound  from  a  marine  propeller  that  is  not  cavitating.  In  particular, 
if  the  shed  vorticity  is  correlated  over  a  substantial  length  along  the  trailing 
edge  of  a  propeller  blade,  singing  occurs .  This  discrete  tone  noise  is  an  im¬ 
portant  contributor  to  both  radiated  and  self  noise  when  it  occurs . 

The  application  of  Kirchhoff’s  formula  to  Eq.  (29)  yields  a  general 
expression  for  radiated  density  fluctuations 
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which  includes  both  Lighthili’s  and  Curie's  theories.  Here  r  =  I  x-yj,  v  is  the 
volume  of  turbulence,  S  is  the  body  surface,  nj  is  a  unit  inward  normal  vector 
to  S  and  the  bracket  [  ]  requires  that  the  quantity  inside  it  be  evaluated  at  the 
retarded  time  t  -  r/c0 . 

If  the  surface  S  is  fixed  and  rigid,  n^Uj  =  0.  Then  the  only  surface 
integral  that  does  not  vanish  is 


d 

dx 


f  r  [Pi]  dS 
Jo 


where  P-  =  njT..  is  the  stress  vector  exerted  on  the  fluid  by  the  surface  boundary 
element  dS .  Tnis  result  forms  the  basis  for  the  determination  of  radiation  from 
finite  rigid  bodies  with  turbulent  wakes,  e.g.,  aeolian  tones. 

Equation  (31)  applies  also  to  the  radiation  from  a  turbulent  boundary 
layer  on  a  rigid  wall.  However,  if  the  wai!  is  of  infinite  extent,  a  result  is  ob¬ 
tained  which  has  been  the  subject  of  controversy  for  ten  years .  The  simplest 
derivation  of  this  result  is  by  means  of  a  reflection  principle  introduced  by 
Powell  (1960).  One  first  notes  that  the  right-hand  side  of  Eq.  (31)  vanishes  if 
the  point  x*  is  interior  to  the  surface  S .  If  S  is  the  infinite  plane  X2  =  0,  it  is 
possible  to  combine  two  expressions  of  the  form  of  Eq.  (31),  one  vanishing  for 
xs  >  0  and  one  vanishing  for  X2  <  0.  This  is  done  by  reflecting  a  turbulent  bound¬ 
ary  layer  flow  in  the  region  X2  >  0  to  form  an  image  flow  in  the  region  xa  <  0. 

It  then  follows  by  symmetry  that  only  a  volume  integral  of  Ty  for  the  real  and 
image  flows  remains,  provided  the  flow  is  assumed  inviscid.  But  this  implies 
that  the  radiation  must  be  of  quadrupole  nature  only,  and  hence  the  intensity 
must  obey  the  law  of  the  eighth  power  of  the  Mach  number . 

An  alternative  proof  of  the  vanishing  of  dipole  radiation  has  been  given 
by  Phillips  (1956b).  A  closely  related  result  concerning  pressure  fluctuations  on 
the  wall  was  obtained  by  Phillips  (1955)  and  by  Kraichnan  (1956)  under  the  assump¬ 
tion  that  compressibility  could  be  neglected  within  the  boundary  layer.  Assuming 
the  turbulence  homogeneous  in  planes  parallel  to  the  wall,  Phillips’  expression  is 
that 


p(x,t)p(x-fr,t)dA(r) 


=  0 
Xj  =  0 

r„  =  0 


(32) 
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where  the  integral  is  taken  over  the  entire  wall  area .  We  may  state  equivalently 
that  the  wavenumber  pressure  spectral  density  vanishes  at  zero  wavenumber  or 
that  the  pressure  correlation  area  on  the  wall  must  vanish . 

Both  of  the  above  related  results  have  important  bearings  on  the  gen¬ 
eration  of  flow  noise  by  turbulent  boundary'  layer  pressure  pulsations .  The  first 
tells  us  that  direct  radiation  from  the  boundary  layer  adjacent  to  a  sonar  dome 
cannot  be  important  and  we  must  look  to  vibration  of  the  dome  itself  for  the 
mechanism  of  flow  noise.  The  second  tells  us  that  there  are  important  restric¬ 
tions  on  the  nature  of  the  forcing  function  causing  the  dome  vibration.  In  par¬ 
ticular,  as  we  shall  see  later,  the  vanishing  of  the  pressure  correlation  area 
weakens  the  validity  of  otherwise  useful  analytical  attacks  on  the  problem  of  a 
plate  or  membrane  subjected  to  random  exciting  pressures . 
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V.  ANALYSES  OF  WALL  PRESSURE  FLUCTUATIONS 


The  fluctuating  pressure  on  a  wall  adjacent  to  a  turbulent  boundary 
layer  is  often  termed  pseudo-sound  in  recognition  of  the  fact  that  it  is  essentially 
non-acoustic  in  nature.  We  have  already  indicated  that  it  is  not  associated  with 
any  significant  far-field  radiation  if  the  wall  is  rigid,  flat,  and  infinite.  How¬ 
ever,  this  pseudo-sound  is  quite  real .  It  is  the  flow  noise  measured  by  a  trans¬ 
ducer  placed  flush  to  a  wall .  It  is  also  the  random  forcing  function  that  sets  a 
sonar  dome  into  vibration  with  consequent  acoustic  radiation  into  the  interior  of 
the  dome.  We  must  therefore  characterize  this  pseudo -sound  in  a  suitable  man¬ 
ner  at  the  outset  of  any  serious  investigation  of  sonar  self  noise. 

The  assumption  that  the  flow  in  the  boundary  layer  is  incompressible 
leads  to  a  great  simplification  in  the  analysis  of  the  relation  of  pressure  fluctua¬ 
tions  to  velocity  fluctuations .  Equation  (29)  then  reduces  to 


dx.  Sx. 

i  1 


d2u.  u. 

P  _ L_i 

°  Sx.  ox. 
»  J 


(33) 


where  PQ  is  the  mean  fluid  density .  By  inspecting  orders  of  magnitude  of  velocity 
and  pressure  terms  in  the  equations  of  motion  taken  very  near  a  flat,  rigid  wall, 
Kraichnan  (1956)  showed  that  the  normal  derivative  of  p  must  be  approximately 
zero  at  the  wall;  i.e.. 


0 


(34) 


If  the  mean  flow  Mach  number  is  small  and  the  wavelengths  under  con¬ 
sideration  are  sufficiently  large,  the  neglect  of  compressibility  within  the  bound¬ 
ary  layer  is  a  reasonable  physical  assumption.  It  means  that  the  purely  hydro- 
dynamic  pressure  fluctuations  dominate  over  the  radiated,  or  acoustic,  pressure 
fluctuations  within  the  iayer.  Not  all  matters  of  interest  can  be  dealt  with  within 
this  assumption,  however.  Ffowcs  Williams  (1964)  pointed  out  that  the  inclusion 
of  compressibility  in  the  analysis  alters  the  conclusion  that  the  pressure  correla¬ 
tion  area  vanishes .  He  found  that  the  correlation  area  is  proportional  to  the 
square  of  the  mean  flow  Mach  number.  For  sonar  applications,  however,  this 
result  only  serves  to  reinforce  the  validity  of  the  incompressible  flow  assumption. 


17 


Using  the  reflection  argument  cited  above,  Kre.ichnan  showed  that 
Eqs .  (33)  and  (34)  imply 


P(x,t) 


b  f  u. 

/  _ Li 

!n  L  >  o  *i 


P 

2n 


dV(y) 

r 


(35) 


where  x  is  a  point  on  the  wall  x2  =  0.  He  also  showed  that  if  the  normal  gradient 
of  the  mean  velocity  is  dui/dy2  then 


By  Sy. 
i  J 


du,  5u'0 

i  2 
^3 


(36) 


where  u'2  is  the  fluctuating  component  of  velocity  normal  to  the  wall .  Hodgson 
(1962)  has  verified  experimentally  that  this  is  the  dominant  term  of  the  velocity 
derivative  products . 

To  proceed  further  with  the  characterization  of  the  statistics  of  wall 
pressure  fluctuations,  we  must  make  certain  broad  conceptual  simplifications 
concerning  the  turbulent  boundary  layer  and  then  test  the  results  of  these  con¬ 
cepts  both  analytically  and  experimentally.  We  consider  coordinates  xx  and  x3 
fixed  in  die  plane  of  a  rigid,  infinite,  flat  wall .  A  turbulent  boundary  layer  is 
established  with  a  uniform  velocity  U  in  the  xx  direction  outside  the  boundary 
layer.  For  a  fixed  position  xx ,  x2  -  0,  x3,  the  ensemble  average 


p(xJtx3;t)  p(x1,x3;  t  +  t) 

is  considered  independent  of  t .  Hence  by  the  ergodic  hypothesis  it  is  equal  to 
a  measurable  time  correlation. 

If  we  can  consider  an  essentially  frozen  pattern  of  turbulent  pressure 
fluctuations  being  convected  along  the  wall  at  some  velocity  Uc  in  the  xl  direc¬ 
tion,  then  with  respect  to  a  moving  coordinate  frame  ST  =  xx  -  Uct,  x2  -  x2 , 
x  =  x,  we  have 

3  ^ 


p(x1?x3;  t)  =  p(x^,x3;  0)  =  p(xx 


U  t,  x  ; 
c  3 


0) 


(37) 
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p(x,,x3:t)  p(xx,x3;t-f  T)  =  p(xlt%;0)  p(x\+ rx  ,x3;0)  (38) 

where  rx  =  -  Uc~ .  For  a  well  established  turbulent  boundary  layer  the  further 
supposition  is  made  that  the  pressure  correlation  is  independent  of  coordinates 
and  x3  . 


We  have  applied  Taylor's  hypothesis,  this  time  to  pressure  fluctua¬ 
tions  rather  than  to  velocity  fluctuations  .  Under  the  above  assumptions,  the 
longitudinal  pressure  correlation  at  the  wail  for  fixed  3c3  =  x3  is 


R(rx  -  U  t,0;0)  =  p(xx  ;0)  p(xx+  r.  -  U  T;0)  =  p(xx  ;t)  p(xx  +  ;t  +  t) 

c  c 


(39) 


It  is  clear  that  the  normalized  form  of  R  will  equal  one  if  t  =  r,  /U  . 

No  analytical  justification  exists  for  the  application  of  Taylor's  hy¬ 
pothesis  to  pressure  correlations.  Moreover,  short  of  experiment,  there  is 
no  physical  basis  for  estimating  the  convection  velocity  Uc  other  than  that 
Uc  <11,  since  a  pressure  correlation  is  the  integral  ef  velocity  correlations  by 
virtue  of  Eq.  (35).  However,  if  we  assume  that  the  wall  pressure  is  statistically 
homogeneous  in  the  moving  coordinate  frame,  it  follows  immediately  that 


p(xj ,  x3  ;t)  p(x j  +  rx ,  x3  +  r3  ;t  +  t)  =  R(rx  -  Uj ,  r3  ;T )  (40) 

j  where  R  is  the  moving-axis  pressure  correlation . 

If  then  it  is  assumed  that 

r 

i 

.  Rfrl  '  =  Rl(ri  '  UcT>  R3(r3>  RM(7>  (41> 

J 

the  measurement  of  the  longitudinal  cross -correlation  of  pressure  (Eq.  40)  with 
|  r3=  0  will  produce  peak  values  of  R  whenever  r,  =  Uc" .  Hence  Uc  can  be  deter  - 

*  mined.  The  envelope  of  these  peaks  is  clearly  the  moving-axis  temporal  auto¬ 

correlation  function  Rm(t)>  which  measures  the  effect  of  decay  of  turbulence  on 
,  the  pressure  correlation.  Taylor's  hypothesis,  of  course,  assumes  that  this 

4  function  is  essentially  constant. 

i 

1 
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The  framework  constructed  above  has  been  used  as  the  basis 
for  several  significant  analyses  of  pseudo-sound.  We  shall  refer  especially  to 
the  studies  of  Kraichnan  (1956),  Lilley  (i960),  and  Ffowcs  Williams  and  Lyon 
(1963).  These  papers  are  strongly  influenced  by  the  fact  that  essentially  all 
earlier  analytical  and  experimental  results  concerning  turbulent  flows  related 
to  velocity  statistics  rather  than  to  pressure  statistics .  Equation  (35),  or  the 
basic  equations  of  motion  from  which  it  is  derived,  form  a  means  for  relating 
pressure  spectra  to  integrals  involving  velocity  spectra .  Both  wavenumber  and 
frequency  spectral  densities  can  be  formed  from  Fourier  transforms  similar 
to  Eq.  (22)  by  use  of  limiting  operations  similar  to  Eq.  (23).  In  parallel  fashion, 
pressure  and  velocity  correlations  can  be  formed,  since  they  are  Fourier  trans¬ 
forms  of  the  corresponding  spectral  densities .  These  analyses  aim  at  reducing 
the  pres  sure -velocity  relations  to  tractable  form  from  which  physically  signifi¬ 
cant  properties  of  the  pressure  field  could  be  inferred.  For  example,  a  spatial 
pressure  correlation  formed  from  Eq.  (35)  leads  immediately  to  two  volume 
integrals  of  fourth  derivatives  of  fourth -order  velocity  correlations.  Clearly 
the  simplifying  approximation  of  Eq.  (36)  is  a  great  aid  to  achieving  any  progress. 

One  might  ask  if  such  effort  is  really  worthwhile.  Why  not  simply 
measure  the  pertinent  properties  of  the  wall  pressure  fluctuations  directly  and 
dispense  with  velocity  considerations?  For,  after  all,  an  adequate  experimental 
description  of  the  wall  pressure  field  alone  would  give  us  the  "forcing  function" 
needed  to  determine  the  vibration  of  a  sonar  dome  and  the  resulting  acoustic 
radiation.  The  answer  to  this  question  lies  in  the  following  considerations: 

a)  The  correlation  area  for  wall  pressure  fluctuations  vanishes  as 
the  wall  b  comes  of  infinite  extent.  Equivalently,  the  wavenumber  spectral 
density  of  pressure  is  zero  at  zero  wavenumber.  There  is  merit,  therefore, 
in  dealing  with  velocity  correlations  with  finite  area  in  this  limiting  case . 

b)  One  can  relate  the  root -mean -square  wall  pressure  to  the  mean 
shear  stress  on  the  wall  by  use  of  the  pressure -velocity  relations.  Since  the 
mean  shear  stress  is  a  measure  of  roughness,  we  have  the  prospect  of  relating 
wall  roughness  to  flow  noise  in  a  systematic  manner. 

c)  The  pressure  correlation  as  expressed  by  Eq.  (41)  is  difficult 

both  to  measure  and  to  interpret .  As  we  shall  see  later  in  our  discussion  of 
experimental  results,  measured  values  of  root -mean -square  pressure  are 
severely  influenced  by  transducer  size,  and  there  is  ambiguity  in  interpreting 
the  relation  of  convection  velocity  Uc  to  the  function  which  measures 

the  decay  of  turbulence.  Recourse  to  pres  sure -velocity  interrelations  is  al¬ 
most  certainly  necessary  if  a  true  understanding  of  the  phenomenon  of  flow 
noise  is  to  be  achieved . 
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(1956)  used  Eq.  (35),  extended  to  include  the  reflected  velocity  field,  to  estimate 
the  gross  properties  of  the  fluctuating  wall  pressures .  He  found  the  following 
proportionality  relations: 


Mean  square  pressure  ~  ~  U4 

Pressure  correlation  length  scale  ~  L  ~  (42) 

Characteristic  frequency  scale  ~  L_1~  U 


where  L  is  the  distance  from  the  point  of  transition  from  laminar  to  turbulent 
boundary  layer  flow  with  free -stream  velocity  U  fixed,  and  the  velocity  propor¬ 
tionalities  are  for  fixed  position  L .  The  apparent  contradiction  implied  by  this 
pressure  correlation  length  relation  with  previously  cited  findings  concerning 
the  correlation  area  is  perhaps  resolved  when  one  considers  that  spatial  homo¬ 
geneity  in  planes  parallel  to  the  wall  is  certainly  not  valid  near  transition.  The 
lack  of  dependence  upon  U  is  a  consequence  of  the  incompressibility  assumption. 
The  characteristic  frequency  scale  expresses  the  dependence  of  the  form  of  the 
pressure  spectral  density. 

Kraichnar.  took  the  frequency -wavenumber  Fourier  transiorm  of  pres¬ 
sure  and  the  right-hand  term  of  Eq.  (33)  and  solved  the  resulting  ordinary  dif¬ 
ferential  equation.  (In  this  and  other  approaches  the  spatial  coordinates  which 
are  usually  transformed  are  those  parallel  to  the  wall,  because  spatial  homo¬ 
geneity  cannot  be  assumed  to  hold  in  the  normal  direction . )  He  was  then  able 
to  obtain  an  expression  for  the  wavenumber -frequency  pressure  spectral  density 
in  terms  of  an  integral  of  the  wavenumber-frequency  spectral  density  of  the 
rigfct-hand  term  of  Eq.  (33): 


*  (k,w) 


2_ 

k2 


e"k(x?  +  x2)  S(x2,x;,k,ui)dx2d4 


(43) 


Here  k=  (k1,k3),  k=jk|,  and  S  is  the  wavenumber-frequency  spectral  density 
of 


-,2 

o  u.u. 

,,  i  J 

o  Sx.  dx. 
i  J 
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pressure  at  the  wall  averaged  over  the  wall . 


k  and  &  then  is  the  mean  square 


Equation  (43)  seems  to  become  singular  as  k-»0  in  apparent  contradic¬ 
tion  to  Phillips’  result,  Eq.  (32).  However,  Kraichnan  and  Ffowcs  Williams  and 
Lyon  demonstrated  that  the  effect  of  the  double  derivative  of  upij  is  to  cause 


£(k,u))--»0  as  k-*0 


The  demonstrations  depend  upon  the  rapidity  with  which  either  velocity  correla¬ 
tions  or  spectral  densities  vanish  as  one  moves  normally  away  from  the  wall. 
Because  of  uncertainty  as  to  the  asymptotic  relations  involved,  none  have  been 
considered  as  complete  proofs.  Recent  experimental  data,  Hodgson  (1962), 
tends  to  support  the  basic  validity  of  the  result . 

Kraichnan  applied  the  simplifying  approximation,  Eq.  (36),  to  Eq.  (35) 
to  obtain  an  estimate  of  the  relation  between  the  root -mean -square  wall  pressure 
and  the  wall  shear  stress: 


Lilley  improved  upon  this  estimate  by  a  somewhat  more  straightforward  analysis 
which  avoided  a  difficult  integration  in  Kraichnan ’s  work  and  obtained 


Tne  essential  assumptions  involved  the  form  of  the  mean  velocity  gradient 
ou.  /ex  and  the  form  of  the  velocity  correlation  R  fr.  ,x_,r_).  Lj.ll-y  took 

X  2  t  22  *  2  3 

the  former  from  Coles  (1956)  experimentally  derived  law  of  the  wake.  He 
chose  the  latter  to  have  the  form  it  would  have  in  isotropic  turbulence.  This 
permitted  the  use  of  Laufer’s  (1954)  measurements  of  the  root-mean-square 
normal  velocity  component  and  Grant’s  (1958)  measurements  of  the  ratio  of  the 
velocity  scale  of  turbulence  to  the  displacement  thickness  of  the  boundary  layer 
to  obtain  the  above  estimate . 
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VI.  MEASUREMENTS  OF  WALL  PRESSURE  FLUCTUATIONS 


The  first  extensive  measurements  of  wall  pre  ssure  fluctuations  were 
those  of  Harrison  (1958).  Small,  flush-mounted  pressure  transducers  were  in¬ 
stalled  in  the  wall  of  a  subsonic  wind  tunnel.  Root -mean-square  pressure,  pres¬ 
sure  spectral  density,  a;.d  longitudinal  and  lateral  pressure  correlations  (using 
pairs  of  transducers)  were  obtained  for  a  tunnel  wind  speed  U  of  from  50  to  200 
feet  per  second. 

Harrison  found  the  root-mean- square  pressure  to  be 

^P2  **  9.5  x  10’  q,  q  =  U2 

2  O 

and  the  convection  velocity  (see  Eq.  (41))  to  be 

U  «  0.8  U  . 
c 

His  measurements  of  pressure  sptictral  density  are  plotted  in  dimensionless 
form  in  Figure  1.  The  experimeroal  form  of  the  pressure  spectral  density 
5(f)  used  by  Harrison  is  related  to  $('u)  of  Eq.  (26)  by  the  relations 

u  =  2^rf  ,  §p(f)  =  4n  $(>n)  . 

Equation  (11)  shows  that  the  displacement  boundary  layer  thickness  6*  is  roughly 
proportional  to  the  distance  from  the  point  of  transition  to  turbulent  flow .  Hence 
Kraichnan’s  proporti.orality  relations,  Eq.  (42),  imply  that 


where  the  form  of  the  function  F  is  given  in  Figure  1.  Harrison's  data  have 
been  replotted  to  confoi  m  to  this  notation. 

Harrison’s  measured  pressure  spectral  density  falls  off  abruptly  for 
t  6*/U  >  1.0.  This  effect  is  due,  in  part,  to  the  failure  of  his  transducer  to 
resolve  the  high-frequency  portion  of  the  spectrum .  The  problem  of  correcting 
for  the  effect  of  finite  transducer  size  is  the  subject  of  several  other  investiga¬ 
tions  which  we  shall  discuss  later. 
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FIG.  I  DIMENSIONLESS  FREQUENCY  SPECTRAL  DENSITY  OF  WALL  PRESSURE 


Harrison's  longitudinal  pressure  correlations  were  obtained  in  nar¬ 
row  frequency  bands.  Thus,  he  actually  measured  the  real  and  imaginary  parts 
cf  the  longitudinal  cross  spectral  density;  i.e.. 


§(ri,0,f)  =  U(rltf)  +  iV(rltf) 


(44) 


where 


U(ri,f)« 


p0<xi.t>  P0<xt  +  ri  jt> 

zi 


V(ri,f)« 


Po(Xl,t)  Po(Xl  +ri’  1  "IT 
_ o 

Af 


p0  signifies  the  filtered  output  signal,  and  a  time  average  is  taken  experimentally. 
The  ergodic  hypothesis  is  employed  to  equal  time  and  ensemble  averages  under 
the  assumption  that  the  process  is  stationary  in  time . 

For  a  sufficiently  narrow  bandwidth,  the  normalized  form  of  U(rx  ,f) 
should  equal 

tTi 

cos  2tvj—  , 
c 

if  the  turbulence  were  frozen.  In  other  words,  U  is  determinable  from  measure- 

c 

ments  of  ,f)  under  Taylor’s  hypothesis. 

Any  deviation  from  a  cosine  form  for  U(ra  ,f)  should  be  indicative  of 
the  decay  of  turbulence.  Following  Bennett  (1956),  Harrison  assumed  the  rela¬ 
tion 


p(xi+r1,t)  =  a(t)+S(t)+Y(t)  (45) 

where  Y(t)  is  that  part  of  the  downstream  wall  pressure  that  is  uncorrelated 
with  the  upstream  pressure,  i.e.. 


p(x1?t)Y(t)  =  0  , 
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anda(l)  and3(t)  combine  with  p{x1,t)  to  yield  U(r\,f)  and  V(r1,f)  respectively. 
He  then  plotted  (see  Bennett  (1957,  p.  30)), 


UlaV  -  ,  .!£! 

Imf  ho 


(46) 


against  Tjf/U  as  a  measure  of  the  coherence  of  p(Xj  +  r^tjwith  p(Xj  ,t).  In 
this  expression,  9^(f)  is  the  frequency  spectrum  of  Y(t). 


Clearly,  the  frequency  transform  of  Eq.  (44)  is  the  longitudinal  cross-correla¬ 
tion  of  pressure,  Eq.  (40),  with  r3  =  0.  The  moving- axis  temporal  autocorrela¬ 
tion  function  Rj^(T )  and  the  measure  of  pressure  coherence 

yo 

i  — i — 


both  describe  the  decay  of  turbulence  in  terms  of  wall  pressure  but  not  in  an 
equivalent  fashion . 

We  have  written  Harrison’s  expressions  in  a  manner  that  would  imply 
that  spatial  homogeneity  is  assumed .  This  was  done  only  to  achieve  conformity 
with  the  approaches  of  other  investigators.  In  point  of  fact,  Harrison  did  not 
assume  spatial  homogeneity  in  the  plane  of  the  wall;  hence  his  approach  is  some- 
whate  more  general  than  those  of  later  investigators . 

Franz  measured  wall  pressure  frequency  spectral  densities  and  rms 
pressures  on  the  submarine  USS  Albacore  over  a  range  of  ship  speeds  and  at 
various  locations  from  near  the  bow  back  to  about  one-quarter  of  the  ship’s 
length  from  the  bow.  His  results  are  reported  in  the  review  by  Richards,  Bull 
and  Willis  ( i960).  They  are  in  general  accord  with  the  wind  tunnel  results  shown 
in  Fig.  1.  The  high-frequency  portion  of  the  spectrum  is  reduced  for  transducers 
located  forward  in  the  thin  portion  of  the  boundary'  layer .  This  is  indicative  of 
the  expected  loss  of  transducer  resolution.  Of  particular  interest  is  the  fact  that 
pressure  spectra  measured  on  a  hydrophone  located  within  the  hull  in  the  forward 
free-flooding  portion  of  the  hull  could  also  be  plotted  non-dimensionally  as  in 
Fig.  1.  The  data  collapsed  to  a  single  spectrum  line  for  all  ship  speeds  except 
the  lowest  speed.  His  results  are  considered  excellent  confirmation  of  the  con¬ 
clusion  that  boundary' -layer  pressure  fluctuations  are  the  cause  of  sonar  self  - 
noise  at  all  ship  speeds  above  that  where  machinery  noise  dominates  and  below 
the  speed  of  onset  of  dome  or  propeller  cavitation. 

Willmarth  (1958),  (1959),  (1965),  Willmarth  and  Wooldridge  ( 1962), 
and  Wooldridge  and  Willmarth  (1962)  carried  out  an  extensive  series  of  measure¬ 
ments  of  wall  pressure  fluctuations.  This  work  is  significant  in  several  respects. 


■%! 
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The  effect  of  the  ratio  of  transducer  radius  to  displacement  boundary  layer 
thickness,  R/6*,  upon  the  wall  pressure  spectral  density  is  clearly  established. 
Figure  1  shows  Willmarth's  (1962)  pressure  spectral  density  measured  beneath 
a  thick  turbulent  boundary  layer  with  R/6*  =  0.166.  We  have  also  plotted 
Willmarth's  1959  spectra  taken  with  large  transducers.  It  is  evident  that  the 
high-frequency  portion  of  the  spectrum  is  lost  for  large  values  of  R/6*. 

Wiilmarth  found  for  a  smooth  wall  that 


5.7  x  10 


-3 


q 


2  .o  "i 
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but  for  a  rough  wall  that 


8 .5  x  10 


q 


His  experiments  involving  roughness  were  very  limited,  but  they  do 
point  out  that  a  significant  increase  in  root-mean-square  wall  pressure  occurs 
when  the  wall  is  roughened.  This,  together  with  the  effect  cf  the  ratio  R/6*, 

may  explain  the  wide  variations  in  ■/pT/q  as  reported  by  various  investigators . 

The  effect  of  roughness  upon  the  ratio  /pa/T  is  less  easy  to  understand  and 
deserves  further  investigation,  particularly  since  the  semi-empirical  results  of 
Kraichnan  and  Lilley  indicate  this  ratio  should  be  independent,  of  roughness . 

Wiilmarth  measured  space-time  wall  pressure  correlations  directly, 
rather  than  through  narrow  frequency  bands.  His  data  can  therefore  be  inter¬ 
preted  in  the  form  of  Eq.  (41).  His  results,  together  with  those  of  Bull,  Wilby, 
and  Blackman  (1963),  are  summarized  in  Fig.  2.  Three  typical  correlations 
are  shown  together  with  the  envelope  R^(" )  of  Eq.  (41).  Bull  found  that  the 
envelope  was  lowered  as  the  Reynolds  number  U5  */v  increased.  Some  difficulty 
in  interpreting  the  longitudinal  correlation  data  arises . 

Wiilmarth  finds  the  convection  velocity  Uc  to  be  dependent  on  trans¬ 
ducer  spacing:  asymptotically,  Uc  —  0.56  U  for  rx  /6*  -•  0  and  Uc  approach¬ 
ing  0.83  U  for  large  values  of  r,  /*•  .  A  physical  explanation  for  this  effect  is 
that  the  large  turbulent  eddies  occupy  the  outer  portion  of  the  boundary  layer 


t Wiilmarth  and  Wooldridge's  (1962)  data  on  rms  wall  pressures  and  frequency 
spectra  were  corrected  by  Wiilmarth  (1965).  We  have  incorporated  the  r  erec¬ 
tions  in  this  report. 
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and  hence  have  a  high  convection  velocity.  The  small  eddies  lie  farther  in  and 
travel  more  slowly.  Since  the  small  eddies  decay  more  rapidly,  they  have  little 
influence  on  correlations  taken  with  large  longitudinal  transducer  separations. 

On  the  other  hand,  for  small  transducer  separations,  the  nearby  small  eddies 
have  not  decayed  appreciably  and  have  a  major  influence  on  the  correlations. 

Thus,  one  should  expect  both  Uc  and  Rj^t)  to  bo  wavenumber- dependent. 

The  longitudinal  spatial  correlation  of  wall  pressure  R(rx  ,0,0)  is 
plotted  in  normalized  form  as  a  function  of  r,/6*  in  Fig.  3.  Measurements 
taken  by  Willmarth  and  Wooldridge  (1962),  Bull,  Wilby  and  Blackman  (1963), 
and  Maestrello  (1965a)  can  be  compared.  Only  those  of  Bull  show  the  negative 
correlation  (for  rx  >  46*)  required  by  Phillips'  theorem.  The  agreement  is 
not  good;  moreover,  no  data  exist  for  transducer  separations  rx  0.66*  . 

In  Fig.  4,  the  lateral  spatial  correlation  of  wall  pressure,  R(0, r3,0), 
is  plotted  in  normalized  form  as  a  function  of  r3/6*.  Here  also  we  may  compare 
the  measurements  taken  by  several  investigators .  All  correlations  are  positive . 
The  early  data  of  Harrison  appears  too  high  and  is  probably  influenced  by 
acoustic  tunnel  noise. 

Maestrello  and  Bull  also  measured  spatial  correlations  at  intermediate 
angles  to  the  mean  flew  direction  and  plotted  isocorrelation  curves  R(r,0) 

=  const,  in  the  rx,  r3  plane.  As  may  be  inferred  from  Figs.  3  and  4,  the  curves 
of  Maestrello  are  elongated  in  the  ri  direction  while  those  of  Bull  are  elongated 
in  the  r3  direction  for  small  correlation  values.  For  large  correlation  values 
near  r  =  0  there  is  a  tendency  toward  isotropy .  However,  the  data  at  small 
transducer  separations  are  too  scanty  for  this  to  be  a  firm  conclusion. 

Considerable  work  remains  to  be  done  in  the  measurements  of  spatial 
correlations  to  verify  the  existence  of  negative  longitudinal  correlations,  to 
establish  the  isotropic  character  of  R(r,0)  near  r*  =  0,  and  to  resolve  the  present 
disagreement  as  to  the  orthotropic  nature  of  R(r,0)  for  large  |  r  j  . 

Corcos  (1962),  (1963)^  (1964)  studied  the  difficulties  involved  in  the 
measurement  of  Uc,  §(uj),  and  p8  by  narrowband  analysis  and  related  his  con¬ 
clusions  to  the  measurements  of  Willmarth  and  Wooldridge  and  other  investiga¬ 
tors  as  well  as  to  his  own  measurements  of  wall  pressure  fluctuations  in  fully 
developed  turbulent  pipe  flow.  He  concerned  himself  particularly  with  the 
measurement  of  spectral  density  §(uu)  and  cross  spectral  density  §(r,  w). 

Writing  the  complex  quantity  i(r,  uu)  as 


LONGITUDINAL  SPATIAL  CORRELATION  OF  WALL  PRESSURE 


LATERAL  SPATIAL  CORRELATION  OF  WALL  PRESSURE 


lie  postulated  that  what  is  actually  measured  is 


*o<->  _  Rell»  . 


R  (0,0) 

O 


$  («1  ) 


=  C(  r,w)  cos(a  -wt) 


(47) 


where  R0(T,t)  is  the  output  cross  correlation  from  a  narrowband  filter  centered 
on  frequency  ou.  Thus 


$  (?,u>)  =  5  (u))  C<  r.o;) 
where  C(G,u>)  =  1  anda(0,x)  =  0. 


ict(r,  u») 


(48) 


Both  C  and  a  are  dimensionless  functions  of  r.ir,  a  convection  velocity 
Uc,  a  turbulence  time  decay  constant  9,  and  some  suitably  defined  longitudinal 
and  lateral  correlation  lengths  ^  If  we  exclude  dependence  on  correlation 
lengths  on  the  basis  of  Phillips’  theorem,  then  for  r3  =  0,  by  dimensional  analy¬ 
sis. 


Corcos  chose  a  = 


(r  w  rs  \ 

-jj — ,  lTT~)  ’ 

/ri®  ri  \ 

ft(r  j,  0»a>)  =  gl—  ,  — -jJ 

'  c  c  '  . 

tTim\ 

=  Al-0— 1  ,  yielding 

/V\  .  r,ui 

\Uc)  6  Uc  ' 


$  (r:  ,  0,a>)  =  $(u!)  A 


(49) 


RQ(r1 , 0,  t) 
Ro(0,0) 


(?) 


cos 


(1) 


U 


(rrucT)  • 


(50) 


He  then  determined  Uc  =  r^x'/a  from  the  ratio  of  the  magnitudes  of  the  real  and 
imaginary  parts  of  the  cross  spectral  density  for  fixed  rx  and  a).  The  function 
A  was  then  determined  as  the  amplitude  of  the  measurement.  Correspondingly 
for  r,  =  0,  the  function 


C(0,r3,o>)  =  ) 


(51) 
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was  determined  and  the  assumption  was  made  that 


/r,u)\  /r  uA 

Kt ) W  ' 


C(rx ,  r3  ,a>)  «  A 


We  may  relate  Harrison's  work  to  that  of  Corcos’  by  noting  that 


C )  • 


%<f> 

5  (£) 


Both  investigators  therefore  assumed  that  the  decay  of  turbulence  could  be  rep¬ 
resented  as  a  function  of  the  dimensionless  ratio  r,  tn/U  . 

1  c 

Applying  the  ideas  of  Uberoi  and  Kovasnay  ( 1953)  on  the  mapping  by 
linear  operators  of  a  random  function  of  several  variables  to  the  problem  of 
the  resolution  of  wall  pressure  by  a  transducer  ol  finite  size,  Corcos  found  that 
the  spectral  density?  m(uu)  measured  by  a  transducer  could  be  represented  as 

$  (uu)  =  f  ?(r;u)6(r)  dr  (53) 

m  J 

where  9  (?)  is  a  function  of  the  transducer  geometry  only  and ? (r,  w)  is  given  by 
Eqs .  (48)  and  (52)  above .  Corcos  then  computed  9  (  r)  for  round  and  rectangular 


transducers,  obtained  the  functions  A 


computed  o  ( rj  1 

©■  ft) 


from  the  measurements 


of  Willmarth  and  Wooldridge,  and  then  obtained  the  ratio  ?  (w)/?(u>  )as  a  func- 

(jur  m 

tion  of  — r -  for  a  round  transducer  of  radius  R  by  numerical  integration  of 
"c 

Eq.  (53). 

Corcos  applied  his  results  to  correct  published  measurements  of 
5  (u))v  in  particular  those  of  Willmarth  and  Wooldridge  for  which  an  excellent 
ratio  R/6*  =  0.166  had  been  obtained.  Even  in  this  case  he  found  a  significant 
lack  of  resolution  of  high  frequencies  with  the  integrated  effect  that 


=  0.61 


P 

where  the  subscript  m  again  denotes  measured 


data. 
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Curcos  pointed  out  that  the  difficulty  in  characterization  of  boundary 
layer  noise  can  be  turned  to  advantage  in  sonar  applications.  A  flush-mounted 
transducer  of  sufficiently  large  radius  will  gain  a  noise  attenuation  equal  to  the 
ratio  The  incoming  acoustic  signal  will  be  affected  only  by  the 

directivity  index  of  the  transducer.  He  gave  the  example  of  ship  speed  «  Uc  = 
20  ft /sec,  R  =  1  ft,  f  =  300  cps  for  which  the  noise  attenuation  factor  was 


$  <uu) 
m 

§(uj) 


0.659  U  2 
_ c 

TTR2  u)2 


2.32  x  10'5 


Bull,  Wilby,  and  Blackman  (1963),  at  the  University  of  Southampton, 
also  carried  out  measurements  of  pressure  spectral  densities  and  space-time 
correlations  in  nai  .ow  frequency  bands.  Bull  (1961)  used  the  assumption  of  a 
separable  cross-correlation  function,  Eq.  (41),  of  an  almost  frozen  field  of 
turbulence  with  constant  spatial  scale  to  reach  the  conclusion  that  A(r1,au)  should 
not  only  be  independent  of  a>  but  also  be  approximately  equal  to  the  moving- axis 

temporal  autocorrelation  function  R^  (\)  .  He  recognized  that  measure¬ 
ments  appeared  to  show  the  contrary  but  pointed  out  that  the  use  of  too  wide  a 
filter  bandwidth  or  the  presence  of  background  noise  in  the  measuring  equipment 
could  lead  to  an  apparent  dependence  upon  uj  .  In  his  1963  paper,  however,  this 

/rjl t>  \ 

viewpoint  was  modified.  He  considered  that  Corcos ’  function  A\  — — /should 

'V 

pertain  to  the  high-wavenumber  components  of  the  pressure  field  which  rapidly 

/r  uj\ 

lose  coherence.  On  the  other  hand,  Al-rp— 1  should  be  replaced  by  a  function  of 
rx  \uc  / 

— — alone  for  the  low-wavenumber  components.  This  view  is  in  general  accord 
uc 

with  our  remarks  based  upon  dimensional  analysis . 


In  Fig.  5  Corcos'  longitudinal  similarity  function  is  plotted  together 

with  Bull’s  narrowband  measurements  of  the  same  quantity.  Corcos  converted 

Willmarth  and  Wooldridge’s  ( 1962)  measurements  of  R»  *(t  )  taker,  in  two  fairly 

M  /‘JIT  \ 

wide  frequency  bands  (300  to  700  cps  and  3000  to  5000  cps)  to  aI  j  by  the 
argument  transformation  \  c  / 


UUl* 


c  c 
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-UNCT 1  ON  FOR  LONGITUDINAL 
TRAL  DENSITY  OF  WALL  PRESSURE 


where  uu  denotes  band  center  frequency.  A  surprisingly  good  collapse  of  the  two 
curves  for  Rj^t)  was  achieved  in  this  way*  Bull's  measurements,  taken  through 
one-third  octave  filters,  show  a  strong  dependence  upon  r}  /h*  at  low  frequencies. 


In  Fig.  6  comparisons  of  Corcos*  and  Bull's  determinations  of  the 


lateral  similarity  function  B 


($)can  1 


be  made.  Corcos  converted  Willmarth’s 


measurements  of  R(0,  r3 , 0)  taken  in  the  above  two  frequency  bands  to 
by  the  argument  transformation 


■ft) 


r  trr 

— 2-  — .  — 3. 

\  LI  * 

3  c 

Similarly  to  the  above  longitudinal  case.  Bull's  more  refined  measurements 
show  a  strong  dependence  on  r3  /o*  at  low  frequencies  . 


Willmarth  and  Rcos  ( 1965)  approached  the  problem  of  pressure  resolu¬ 
tion  from  a  different  viewpoint  from  that  cf  Corcos.  They  sought  an  expression 
for  2(w)  in  terms  of  ^m(r,w)  and9("r),  or  transforms  thereof,  and  then  used 
Corcos*  similarity  hypothesis  and  empirical  relations  to  represent  r,a0 
rather  than  §(*r,a.)-  Using  the  wavenumber  transform  *{  k)  of  the  transducer 
geometry  function  SCx) ,  they  obtained  the  relations 


where 


m 


(2tt  y 


f  «k, -)  t(k). 

(54) 

f  <?u,)  dk. 

m  i(k) 

(55) 

f  R(  r,  t)  e  i  k  rdr  . 

A  curious  circumstance  resulted.  For  a  circular  transducer  of  radius  R  , 


where  Jx(kR)  is  the  Bessel  function  of  first  kind  of  order  one.  At  wavenumbers 
corresponding  to^  the  zeros  of  J  the  integrand  of  Eq.  (55)  becomes  infinite  for 
the  form  of  §m(  k,  in)  following  from  Corcos'  assumption  applied  to  §m(*r, 
instead  of  to  $(  r,  a).  However,  the  frequency  transform  of  Eq.  (54)  indicates 
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that  »  (  k,oj)  has  precisely  the  same  zeros  as  ?(  k).  The  difficulty  with  Corco s' 
model  applied  to  *  (  r,x)  is  that  the  measurements  upon  which  the  similarity 
functions  A  and  B  are  based  were  not  taken  with  sufficiently  small  transducer 
separation . 


Willmarth  and  Roos  obtained  a  new  estimate  for  ;m(a’)/r(>Ju)  using  the 
similarity  functions  A  and  B  in  Eq.  (51)  but  truncating  the  integration  short  of  the 
first  zero  of  y(  k).  The  effect  of  this  truncation  was  to  make  *m(u>)/$('x)  dependent 
upon  R/b*  as  well  as  upona>R/Uc.  By  graphical  integration  of  Eq.  (55)  they 
found  *m(w)/*(a)to  correspond  to  Corcos’  results  for  xsR/Uc  5  1.  For 
^R/l5c  >  1  their  results  show  Corcos*  correction  to  be  too  great.  Willmarth 
has  reported  informally  that  he  has  tried  to  apply  Eq.  (55)  and  Corcos*  model  to 
correct  his  own  experimental  data  for  four  different  sizes  of  transducers  (see 
Willmarth  and  Wooldridge  (1962)).  The  corrections  made  the  larger  transducers 
yield  higher  values  of  £(x)at  high  frequencies  than  the  small  transducers. 


A  criterion  for  the  validity  of  correlation  models  was  proposed  by 
Ffowcs  Williams  (1960)  for  velocity  correlations  and  applied  by  Cbandiramani 
(1965)  to  wall  pressure  correlations.  Suppose  R(r,  t)  is  a  one-space  dimensional 
cross -correlatioji  measured  in  a  reference  framejnoving  at  the  convection 
velocity  U  .  Then  for  sufficiently  small  values  r  and  t. 


R(7,-)  ^  K(0,-) 
R(r,0) 


(56) 


Inequalities  (56)  in  fact  uniquely  define  a  convection  velocity  associated  with  the 
broadband  cross-correlation  function.  For  if  a  ridge  of  local  maxima  of  the 
cross -correlation  measured  in  any  reference  frame  exists,  then  it  follows  that 
there  is  only  one  reference  frame  moving  at  Uc  in  which  this  ridge  cannot  be 
detected.  Under  the  assumption  that  the  wall  pressure  is  a  stationary  random 
function  of  space  and  time  in  any  reference  frame,  it  follows  that  if  r  =  r-Uct, 
the  cross-correlation  K  measured  in  a  fixed  reference  frame  is  related  to  R  by 

R(r,  ’)  =  R(r  -  U  "  )  .  (57) 

c 

It  then  follows  that  for  sufficiently  small  r, 

R(0,")  *  R(U  -,') 
c 

*  R(U  -,0)  (58) 

c 
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Chandiramani  (1965)  pointed  out  that  certain  models  of  the  cross- 
spectral  density  used  in  recent  analyses  of  experimental  measurements  of  wall 
pressure  do  not  satisfy  inequalities  (58). 

Bull  (1961)  used  the  approximate  relation  for  r3  =  0, 

}(*i  •«•>  -  rm  (  if 

If  one  takes  the  frequency  transform  of  this  relation,  one  finds  that 


R(U  t,0)  <  R(0,  T) 
c 

in  contradiction  to  FfowcsWilliams*  criterion.  One  may  readily  show  that 
Corcos'  (1962),  (1963),  (1964)  model  of  the  cross -spectral  density  also  fails 
to  satisfy  the  criterion  if  the  criterion  is  assumed  to  hold  for  cross  correlations 
measured  in  narrow  frequency  bands. 


Expressions  (56)  and  (57)  lead  to  the  more  general  inequality 
R(r,  T)  ^  R(Ut,t)  . 

Chandiramani  pointed  out  that  this  inequality  implies  that  the  proper  method  of 
determining  Uc  from  R(r,')  experimentally  is  to  maximize  R  by  varying  r  with 
~  fixed.  Most  determinations  of  Uc,  however,  have  been  mace  by  fixing  r  and 
varying  t.  That  these  two  methods  are  not  in  general  equival  ?nt  is  evident 
from  the  example  of  the  separable  correlation 


R(r,"  )  =  R.(r-U  -)R.,(T) 

-  C  i\l 


foj*  which 


o  R  n  u  .  Sr  _  e>  *  /-\ 

—  =0  but—  =R.,(  ) 

Sr  ..  _  =  '  . .  _  M 

r  =  U  r  =  U 

c  c 


Since  Rj^j  (T)  <  0  for  "  >  0,  the  latter  method  yields  too  large  a  value  for  Uc. 

The  determinations  of  Uc  from  cross -spectral  densities  made  bv  Harrison  and 

Corcos  suffer  from  essentially  this  same  deficiency. A 

5  It  appears  that  the  first  of  the  inequalities  of  Eqs.  (56)  and  (58)  as  well  as  the 
inequality  of  (58a)  are  valid  whenever  a  convection  velocity  Uc  can  be  established, 
'’he  second  and  stronger  inequalities  of  Eqs.  (56)  and  (58)  are  only  plausible, 
lor  it  is  possible  to  construct  space-time  correlations  for  which  these  in¬ 
equalities  do  not  hold. 
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The  error  in  0C  determinations  by  the  usual  experimental  techniques 
is  not  of  itself  particularly  important  for  Rj^(T)  if  the  decay  of  turbulence  is 
sufficiently  slow.  Thus  Corcos'  choice  ofa(rx,u))  =  r1a/Uc  is  not  a  gross  ap¬ 
proximation.  On  the  other  hand,  the  use  of  the  similarity  function  A(r  ,uj)  = 
A^cu/UjJ  leads  to  the  more  fundamer'ial  difficulty  that  the  resulting  corrections 
to  frequency  spectral  densities  measured  with  transducers  of  finite  size  are  too 
large  at  high  frequencies.  Chandiramani  proposed  the  alternative  concept  of  a 
narrowband  cross -correlation  measured  through  a  (hypothetical)  wavenumber 
filter.  One  thus  obtains  instead  of  Eq.  (48)  the  relation 

-ik,U  t 

=  U^)  A(kxUcT)e  *  C 

where  A(ktUcT)  is  essentially  R^jCO  but  for  the  wavenumber  component  kx  only. 
No  difficulties  with  Ffowcs  Williams'  inequalities  arise  from  this  approach. 

One  can  utilize  existing  experimental  data  simply  by  making  appropriate  conver¬ 
sions 


Chandiramani  used  the  following  generalization  of  the  wavenumber -frequency 
transform  of  Eq.  (41)  and  writes  the  spectrum  as 

t(k,a)  --  o(k)  $  ,(k,U  (k);o»-k.U  ) 

Me  1  c 

where  \ .  is  the  wavenumber  transform  of  R. .(").  In  turn, 

M  M 


-k  U  t 
l  c 

R..O  =  A(kiU  ")  =  e 

where  &  =  9  from  Willmarth's  measurements  of  R(rlT  ~ ).  He  assumed  §(k)  to  be 
isotropic  in  k,  but  to  be  related  to  a  finite  correlation  radius  (an  invalid  assump¬ 
tion  by  Phillips'  theorem).  He  then  calculated  the  correction  to  $(uj)  for  finite 
transducer  radius  in  the  manner  of  Corcos  but  found  the  correction  factor  for 
high  frequencies  to  be  much  less  significant  than  Corcos  found  it  to  be.  His 
method  permitted  the  ratio  of  transducer  radius  to  displacement  boundary  layer 
thickness  to  influence  the  correction  in  general  accordance  with  experimental 
trends.  We  note,  however,  that  this  ratio  is  introduced  into  Chandiramani’s 
analysis  in  an  entirely  different  fashion  from  the  manner  in  which  it  appeared  in 
that  of  Willmarth  and  Roos . 

Chandiramani  followed  Corcos  in  that  he  assumed  the  form  of  the 
true  *( k,  ju>,  but  he  did  not  accept  Corcos'  similarity  variables.  Willmarth  and 
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Roos  used  the  similarity  variables  but  applied  them  to  the  form  of  the  measured 
§  (  r,o).  It  would  appear  desirable  to  adapt  Chandiramani's  approach  to  measured 
data.  This  would  require,  however,  an  experimental  realization  of  the  concept 
of  a  narrowband  wavenumber  filter . 

Chandiramani  made  the  interesting  observation  that  a  transducer  of 
finite  area  should  measure  essentially  the  moving- axis  frequency  spectrum 

at  high  frequencies  rather  than  the  fixed-axis  frequency  spectrum  *(<jj). 

An  experimental  confirmation  of  this  analytical  conclusion  would  be  difficult 
because  of  internal  high-frequency  noise  problems  in  measuring  equipment.  It 
would  be  desirable  nevertheless.  For  as  we  shall  see  later,  ^(a>)  or,  equiva¬ 
lently,  the  temporal  decay  of  wall  pressure  turbulence  is  an  important  quantity 
in  determining  structural  response . 

Gilchrist  and  Strawderman  (1965)  measured  the  responses  of  two 
piezoelectric  transducers  in  one -third  octave  frequency  bands  to  boundary  layer 
pressure  fluctuations  in  a  pipe  flow  facility  over  a  range  of  convection  velocities. 
The  effective  radius  of  each  transducer  was  determined  by  measuring  its  voltage 
response  to  the  loading  of  a  weighted  pin  located  at  successively  increasing  radii 
from  the  center.  The  ratio  of  the  spectral  densities  of  the  two  transducers  should 
equal  the  ratio  of  Corcos*  (1963)  correction  factors,  with  all  quantities  expressed 
in  terms  of  iur/U  .  The  ratio  of  effective  radii  was  2.56,  and  the  spectral  den¬ 
sities  were  determined  over  a  range  .01  <  &r/Uc  <  4. 

Gilchrist  and  Strawderman  constructed  an  experimental  correction 
factor  for  transducer  diameter  by  an  iterative  use  of  their  data  for  wr/Uc  from 
.01  to  4.  The  agreement  with  Corcos'  semi-empirical  correction  factor  was 
within  experimental  error  over  this  dimensionless  frequency  range.  It  appears, 
however,  that  this  confirmation  of  Corcos'  result  is  not  valid  for  transducer 
radius  corrections  significantly  greater  than  the  2.56  ratio  used  in  their  experi¬ 
ment.  Corcos'  similarity  hypothesis  is  implicit  in  their  work  as  well  and  hence 
was  not  tested. 

It  would  be  premature  to  assess  the  relative  merits  of  the  several 
viewpoints  currently  held  concerning  the  proper  statistical  representation  of  the 
wall  pressure.  In  spite  of  the  high  degree  of  refinement  of  measurement  tech¬ 
niques,  basic  questions  remain  concerning  the  interpretation  of  experimental 
results.  We  have  emphasized  only  the  situation  of  a  fully  developed  turbulent 
boundary  layer  adjacent  to  a  smooth,  flat  plate.  The  effects  of  surface  rough¬ 
ness,  pressure  gradient,  flow  separation,  and  transition  remain  largely  un¬ 
explored  . 
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VII.  THE  INFLUENCE  OF  STRUCTURE  ON  FLOW  NOISE 


vA 

Except  for  the  special  case  of  flush -mounted  transducers,  all  sonars 
are  mounted  inside  free -flooding  domes.  The  domes  are  usually  constructed  of 
thin  steel  plates  welded  over  a  rib-like  frame.  A  few  domes  have  been  con¬ 
structed  of  fiberglass  or  rubber;  we  shall  not  discuss  these  specifically,  although 
our  remarks  in  general  apply  to  such  domes  as  well. 

There  is  some  point  in  reminding  ourselves  from  time  to  time  as  to 
the  purpose  of  a  sonar  dome.  It  iS  to  permit  sonar  operation  at  speed;  that  is, 
to  eliminate  the  violent  self  noise  that  would  occur  from  cavitation  and  turbulence 
resulting  from  flow  about  the  non -streamlined  sonar  elements.  That  there  is  no 
other  reason  has  been  demonstrated  in  at  least  one  instance  where  low-speed 
sonar  performance  was  actually  improved  by  removing  the  dome . 

Despite  the  purpose  of  the  dome,  its  design  has  heretofore  been  based 
on  criteria  not  directly  associated  with  flow -noise  reduction ,  Its  hydrodynamic 
form  is  dictated  by  drag  and  cavitation  considerations,  and  its  structural  form  is 
a  compromise  between  strength  and  acoustic  signal  transmission  requirements . 
Domes  have  not  been  designed  to  minimize  flow  noise,  because  no  quantitative 
criteria  for  the  effects  of  smoothness,  form,  and  structure  on  internal  acoustic 
radiation  from  boundary  layer  pressure  pulsations  have  been  available . 

Turbulent  boundary  layer  pressure  pulsations  excite  principally  bending 
waves  in  the  sonar  dome  structure.  These  bending  waves,  in  turn,  cause  acoustic 
radiation  both  outside  and  inside  the  dome .  The  internal  radiation  is  detected  at 
the  sonar  transducer  as  self  noise  from  flow.  In  reality,  both  the  vibratory  mo¬ 
tion  of  the  dome  and  the  internal  acoustic  field  are  reverberent;  hence,  the  actual 
self  noise  measured  by  a  transducer  element  is  a  spatial  average  over  all  wave 
fronts  to  which  die  element  is  responsive.  Both  near -field  characteristics  of 
plate  radiation  and  transducer  directivity  must  play  an  ultimate  role  in  a  calcu  ¬ 
lation  of  the  self  noise  received . 

Because  of  the  complexity  of  the  total  problem,  most  analytical  investi¬ 
gations  have  been  restricted  to  the  determination  of  far -field  radiation  from 
either  infinite  or  bounded  flat  plates  excited  by  boundary  layer  turbulence. 

Ribner  (1956)  and  Corcos  and  Liepmann  (1956)  studied  the  case  of  the  infinite 
plate.  Kraichnan  (1957)  gave  the  earliest  treatment  of  the  case  of  a  periodically 
supported  plate.  A  significant  limitation  of  this  work  is  that  the  effect  of  acoustic 
loading  on  plate  vibration  is  largely  neglected.  Such  an  approximation  is  valid  in 
air,  where  mass  loading  of  the  air  on  the  plate  is  negligible  and  the  damping  from 
acoustic  radiation  is  small  compared  with  the  internal  or  applied  structural  damp¬ 
ing  of  the  plate.  In  water,  however,  neither  of  these  effects  is  necessarily  small, 
and  the  problem  is  thus  much  more  difficult. 
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The  problem  for  an  infinite  plate  in  water  can  be  formulated  as 
follows:  The  equation  for  small  normal  displacement  r|  (xa  ,  x3  ,  t)  of  a  plate 
excited  by  wall  pressure  p  (xx  ,  x3  ,  t)  is 


D  '74n  +  p  h  r| 
p  tt 


p_  “Vp 


(59) 


where  D  = 


12 


Eh3 

(1-os) 


is  the  flexural  rigidity  of  the  plate,  h  is  the  plate 


thickness,  Pp  is  the  plate  density,  and  E  and  o  are  Young's  modulus  and  Poisson’s 
ratio,  respectively.  Structural  (hysteresis)  damping  of  the  plate  can  be  crudely 
accounted  for  by  writing  E  (1  +i  s)  in  place  of  E,  where  e  is  the  material  loss 
factor .  Alternatively,  one  can  insert  an  arbitrary  Rayleigh  damping  term 
Pph  3"nt  in  the  left  hand  side  of  Eq.  (59). 


The  terms  P+  are  the  acoustic  radiation  pressures  above  and  below  the 
plate.  The  radiation  field  is  related  to  the  plate  vibration  by  the  boundary  con¬ 
dition  . 


and  die  Bernoulli  equation 


a_n 

at 


a  0 


p  =  -  p 


a  0 

o  at 


(60) 


(61) 


where  PQ  is  the  fluid  density  and  0(x1 ,  x2  ,  X3  ,  t)  is  the  fluid  velocity  potential. 
The  fluid  velocity  potential  must  satisfy  the  wave  equation 


?20  -  -4-  5>  =  0  (62) 

c2  tt 
0 

where  cQ  is  the  speed  of  sound  in  the  fluid . 

By  taking  Fourier  transforms  of  the  above  relations  we  can  write  the 
wavenumber -frequency  spectral  density  of  the  field  pressure  P  at  a  distance 
Xa  from  the  plate  as 


(63) 
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where y=  %'k2  -k*  ,  kQ  =  x/cq  »  *  (k*,  oj)  is  the  wavenumber -frequency  spectral 
density  of  the  wall  pressure,  and  k^  =  where  the  plate  wave  speed  Cp 

satisfies  the  dispersive  relation 


By  integrating  Eq.  (63)  over  k  .  the  frequency  spectral  density  of  field  pressure 
may  be  obtained.  A  further  integration  over  x  will  then  yield  the  mean  square 
field  pressure.  Both  quantities  derived  in  this  manner  are  averages  over  the 
plate  surface.  Thus  the  problem  of  predicting  the  radiation  field  at  a  distance 
x2  from  a  plate  undergoing  turbulent  boundary  layer  excitation  described  by 
i  (k,  x)  can  in  principle  be  solved  explicitly.  The  effect  of  fluid  loading,  indi¬ 
cated  by  the  term  v  ,  adds  a  severe  complexity,  and  no  general  solutions  have 
as  yet  been  obtained . 

We  may  note  from  Eq.  (63)  that  when  k  <  kQ,  Re'*  =  0  and  hence  no 
attenuation  of  with  increasing  x2  occurs.  It  is  also  evident  that  when  k  -  kp 
we  should  expect  a  marked  increase  in  the  magnitude  of  £ f.  Since  k0  is  usually 
significantly  less  than  kp  at  frequencies  x  of  interest  in  sonar  applications,  the 
previously  cited  result  that 

?  (k  .  x)  — - — ►  0  as  k  — -  0 

implies  that  the  dominant  contribution  to  the  radiation  field  comes  from  wave- 

numbers  near  k  . 

P 

A  remark  is  in  order  concerning  the  representation  of  the  wall  pres¬ 
sure  spectrum  t(lc ,  x)  under  the  assumption  of  convected  frozen  turbulence. 
From  Eq.  (41)  with  R-  ( ’)  =0  (no  decay)  the  pressure  correlation  can  be 
written  as 


p  (x: ,  x3 ;t) p  (x,  +r-; ,  x3  +  r3 ;  t  +  r )  =  R  (  r;  '• )  =  R^  (r:  -  U  " ,  r3  )  .  (65) 

Note  that  it  is  not  necessary  to  factor  R0  into  R.  ( r.  -  U  t)  R3  (r3  ). 

We  have  the  basic  relation 


•  \ 
/ 


-  ik  •  r 
e 


d  * 


(66) 
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Introducing  Eq.  (65)  into  (66)  and  changing  variables  of  integration  to 
~  ri  "  UcT,  =  T,  we  can  show  directly  that 


Mk,  *)  =  5(‘JJ'Uck x3;0)  e'1  (kl  ^+k5X3>  d§  dx3 
where  6  is  the  Dirac  delta  function.  But 


/ 


4>(k,  a)  e  iv  dai  =  -  J  R(r,  T)e  ik  rdr 


hence,  setting  t  =  0  in  the  relation  we  obtain 


?  (k,  a)  =  5  (a  -  UJq  )  j  ?  (k,  ju)  d 


=  6  (vju  -  U  kj)  ?  (k) 


(67) 


Eq.  (67),  or  modifications  of  it  for  known  R  ^  (t)  f  1,  form  a  basis 
for  the  spectral  representation  of  the  turbulent  wall  pressure  as  a  forcing  function 
on  the  plate . 


Ribner  (1956)  dealt  with  the  radiation  from  an  infinite  flat  plate  sub 
jected  to  a  convecting  frozen  field  of  turbulence.  He  ignored  the  effect  of 
radiation  loading  upon  plate  vibration;  hence,  instead  of  Eq.  (63)  he  obtained 
the  equivalent  of  the  relation 


»£  (X2  ,  k,  a.') 


-2Re  yxa  .  _ 

e _ i  (k,  a) 


/  p  }l\2  V/ 

L  .4  1 

M  - 1 

v  po '  IV 

V  J 

\  i 


(68) 


For  wall  pressure  convecting  at  velocity  Uc  without  decay 
R  (r,  * )  =  Rq  (r:  -  U  t  ,  r3) 

for  which 


J(M)=  5(k,U  -x)  ?  (k) 


by  Eq .  (67) .  Then 
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(69) 


$f(x2. 


‘f 


(x2,  k,  uu)  dk 


A  far-field  pressure  frequency  spectrum,  lim  §£  (x2,  w),  is  finite 

x2  -  ® 
uu 

only  if  Rey  =  0.  However,  if  Uc  <  cq  and  kx  =  —  ,  this  cannot  occur  for  any  k3  . 

Thus  an  infinite  plate  subjected  to  a  subsonic  frozen  convected  pattern  of  wall 
pressure  cannot  radiate  to  the  far  field . 


Ribner  therefore  chose  to  compute  j* $£  (0,  k)  dk  ,  the  mean  square 

field  pressure  at  the  plate  surface  averaged  over  the  plate  surface.  He  intended 
to  assume  isotropy  of  the  wall  pressure  in  the  plane  of  the  wall  as  expressed  by 

R  (r,  0)  =  r,  r  =  j  r  j  . 

Here  J*R(r,  0)  dr  =  p2  2ttX2,  where  X  is  a  correlation  radius.  However,  he 
inadvertently  used  a  three-dimensional  transform 


*(k)  =  2 


kX3 

(1  fk2X2)2 


(70) 


instead  of  the  proper  two-dimensional  transform 


•>2 

•\ 


(1  +kaXs)3/  2 


•(71) 


as  his  wavenumber  spectrum.  He  then  introduced  structural  hysteretic  damping 
by  replacing  kp  by  kp  (1+s),  where  e  is  the  material  loss  factor.  The  effect 
of  the  plate,  then,  is  to  filter  $  (k)  in  a  narrow  wavenumber  band  about  k  =  kp. 
His  calculations,  based  on  the  assumptions  that  X  is  proportional  to  6*  and 
Uc  %  U,  led  him  to  the  conclusions  that  for  x2  =  0, 
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P2  =f  $f  (0,  k)  dk  ~ 


u5  6*3 

for 

U  6* 

h° 

C .  h 

U3  6* 
h 

for 

U  6* 

cih 

(72) 


large  . 


Here  = 


E 


P  (1  -o2) 
P 


1 

J 


is  the  longitudinal  (bulk)  speed  of  sound  in  the  plate . 


In  light  of  later  investigations,  Ribner’s  presumably  accidental  choice 
for  §  (k)  was  fortuitous.  For  as  given  by  Eq.  (70),  I  (k)  —  0  as  k  —  0  in 
accordance  with  the  requirements  of  Phillips'  theorem .  Had  Ribner  used  the 
$(k)  given  by  Eq.  (71)  corresponding  to  rigidly  convected  isotropic  wall  pressure, 
he  would  have  found 


P2 


U4  6*a 
ha 


as 


U  6* 
C^h 


0  , 


U6* 

C^h 


(73) 


Corcos  and  Li epmann  (1956) also  treated  the  case  of  the  infinite  plate. 
They  used  Rayleigh’s  classical  formula 


P  (x,  t)  = 


(74) 


for  the  determination  of  P2  in  terms  of 

bounded  region  of  plate  vibration  is  being  considered  by  requiring  x2  <  <  R  where 
|y  j  ^  R  is  the  domain  of  integration  in  Eq.  (74).  However,  the  admittance  function 
they  use  is  essentially  the  same  as  that  of  Ribner,  Eq.  (68),  and  is  devoid  of  any 
finite  plate  resonances .  Their  work  differs  from  that  of  Ribner  in  that  they  are 
led  from  Eq.  (74)  to  consider  a  correlation  area  Xs  of  plate  acceleration  given  by 


( 


a2-n\  _ 

d ~r  j  .  They  give  the  impression  that  a 


t(0,  0)  >: 


where 


■■ 


(r,  0)  dr 


d2  q  —  c)2  ti  —  -• 

(r* T)=  n*  iy’l)  ^  (y-rr»t  +  7) 


dt2 
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From  this  they  obtain  the  approximate  relation 


at  a  point  x  many  X-lengths  away  from  the  plate. 


They  next  assume  a  rigidly  convected  pattern  of  isotropic  wall  pres¬ 
sure  with  correlation  radius  proportional  to  6*.  The  "radius"  A.  is  found  to  be 
proportional  to  Cjh/U  by  utilizing  the  fact  that  only  wavenumbers  k  =  kp  con¬ 
tribute  significantly  to  Ps  .  Small  Rayleigh  damping  is  included  in  the  plate 
equation,  and  the  result 


(76) 


is  obtained.  The  function  f  is  not  determined  explicitly,  as  they  did  not  specify 
the  form  of  §(k).  It  is  clear,  however,  that  neither  Ribner’s  relations, 

Eq.  (72),  nor  those  of  Eq.  (73)  can  be  obtained  from  Eq.  (76).  One  suspects 
that  the  approximation  of  Eq.  (75)  is  at  fault. 

Kraichnan’s  (1957)  analysis  differs  from  that  of  Corcos  and  Liepmann 
in  several  respects,  not  the  least  of  which_is  its  complexity.  He  attempted  to 
infer  more  specific  information  about  §  (  k  ,  oj)  from  the  theory  of  boundary 
layer  pressure  pulsations.  He  assumed  that  the  infinite  plate  was  made  up  of 
a  flat  array  of  square,  mechanically  independent  sections  of  side  L  satisfying 
boundary  conditions  of  zero  displacement  and  zero  bending  moment  at  the  edges. 


Two  features  of  Kraichn,  a’s  assumed  wall  pressurejspectrum  are 
worthy  of  note:  First,  the  frequency  spectrum  (integral  of  ~(k  .  uu)  over  k)  is 


assumed  to  cut  off  at 


of  the  order  of  one.  Second,  he  noted  that  the 


approximation  of  Eq.  (36)  neglects  transverse  contributions  to  the  wall  pressure 
spectrum  which  might  be  dominant  at  low  wavenumbers.  The  resultant  form  of 
the  frequency  spectrum  which  he  assumed  does  not  correspond  well  with 
measurements;  hence,  the  conclusions  he  drew  from  his  analysis  are  subject 
to  this  limitation . 


Although  Kraichnan  also  neglected  the  effect  of  fluid  loading  on  plate 
vibration,  he  did  make  a  significant  contribution  by  considering  the  effects  of 
resonant  modes  of  finite  plate  vibration  under  the  assumption  of  a  Rayleigh 
damping  in  each  mode.  Equation  (59)  with  the  right-hand  side  set  equal  to  zero, 
together  with  the  boundary  conditions  of  zero  displacement  and  zero  bending 
moment  at  die  edges  of  a  square  plate  of  side  L,  yields  the  orthonormalized 
eigen  functions 
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for  integer  m  and  n . 


(77) 


The  corresponding  eigenvalues  or  natural  frequencies  uumn  of  the  plate 
are  given  by  the  relation  _ 
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mn 
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We  note  in  Eq.  (78)  that  the  resonant  frequency  u^nn  is  proportional  to 
the  square  cf  the  mode  numbers  m  or  n. 


To  obtain  the  vibratory  response  of  a  finite  plate,  Kraichnan  expanded 
the  temporal  Fourier  transforms  of  the  normal  wall  velocity  and  the  wall  pressure 
in  terms  of  the  eigenfunctions;  i.e.. 


V(X,  ID)  =  -iU)Tl(x,  0!)  =  £  v  (u>)n  (x) 

mn  mn 

m,  n 


m.n 


(79) 


Tlie  modal  admittance  is  determined  from  the  transform  of  the  plate 
equation  as 
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(80) 


From  the  orthogonality  of  the  eigenfunctions  and  the  wall  pressure  cross 
spectral  density  relation 


§(r,  uj)  =  lim  (2n)  p  (x,  uu)  p  (x  +  r,  -  u>)  f  $(i<,  uu)  e*k  r  die 

T— 06  - - -  J  (81) 


the  modal  frequency  spectral  density  of  wall  pressure  can  be  derived  as 


S  (w)  =  lim  (2tt)  F  (w)  F  (■«») 
mn  _  mn  mn 

T  — 35  - 


lotr4  m2  n2 
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(k,  uu)dk  . 
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From  Eq.  (80)  it  follows  that  the  modal  frequency  spectral  density 
of  normal  wall  velocity  is 


V  (^)  =  lim  Vmn(J,>  Vmn  (  -  a<)  = 
nin  _  mn  mn  mn 

1  —  ®  - — - 


mn 


(83) 


Certain  physical  facts  are  evident  from  Eqs .  (82)  and  (83) .  The  modal 
pressure  frequency  spectrum  2mn(ju)  is  clearly  controlled  by  wavenumbers  coin¬ 
cident  with  the  modal  wavenumbers.  It  then  follows  that  the  modal  velocity 
spectral  density  Vmn(ai)  is  peaked  at  the  plate  resonant  frequencies.  This  close 
interrelation  between  wavenumbers  and  frequencies  could  perhaps  have  been 
anticipated  by  an  inspection  of  Eq.  (78),  expressing  the  frequency-wavenumber 
relation  for  free  standing  waves  on  a  finite  plate.  The  cosine  terms  in  Eq.  (82), 
however,  represent  a  "smearing"  of  the  response  over  sones  in  kak3  wavenumber 
space,  each  centered  on  the  discrete  points  representing  the  wavenumbers  of  the 
eigenfunctions.  By  comparison,  the  wavenumber -frequency  relation  for  an  in- 

k*  =w  J. 

P 

the  field  pressure  response  as  shown  in  Eq.  (63). 


finite  plate  is  the  continuous  relation  k 


=  k2  = 


which  controls 


One  cannot  proceed  further  with  the  analysis  of  the  vibratory  response 
of  a  finite  plate  or  with  the  determination  of  the  resultant  field  pressure  wave¬ 
number-frequency  spectral  density  without  making  an  important  additional 
physical  assumption.  It  is  clear  by  Eq.  (79)  that  if  one  were  to  attempt  to  com¬ 
pute,  say,  the  frequency  spectral  density  of  normal  velocity  of  the  plate,  contri¬ 
butions  of  a  very  complex  nature  would  arise  from  coupling  between  different 
plate  modes.  Kraichnan  made  the  assumption  that  contributions  to  the  vibratory 
response  and  to  the  field  pressure  from  different  plate  modes  are  statistically 
uncorrelated.  He  further  assumed  that  contributions  to  the  field  pressure  from 
different  vibrating  plate  elements  are  statistically  uncorrelated. 

One  can  demonstrate  formally  that  the  contributions  of  different  normal 
modes  of  a ^quare  plate  to  the  wavenumber -frequency  normal  velocity  spectral 
density  V  (k,  jj)  are  statistically  uncorrelated  in  the  special  case  where  the  wall 
pressure  correlation  function  R(r,  ")  is  written  (in  contradiction  to  Phillips’ 
theorem!)  as 

R(r,  ~)  =  p2  A  5(ri)  5  (r3)R.(')  (84) 

where  A  is  the  correlation  area.  The  cross  spectral  density  of  wall  pressure 
corresponding  to  this  correlation  also  contains  the  Dirac  delta  functions 
5(r:)  &(r3).  It  is  this  fact  that  permits  the  use  of  the  orthogonality  of  the  normal 
modes  in  the  derivation  of  the  result 
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(85) 


Clearly,  each  mode  contributes  independently  to  the  velocity  spectrum . 
The  modal  admittance  |Ymn(tn)  |2  is  obtained  from  Eq.  (80).  The  modal  fre¬ 
quency  spectral  density  of  wall  pressure  Jmn(ti))  is  given  by  Eq.  (82).  For  the 
particular  correlation  function  chosen  in  Eq.  (84)  it  can  be  shown  that  §mn(u>)  is 
independent  of  mode  numbers  m  and  n. 

In  Eq.  (84)  the  convection  velocity  Uc  is  assumed  zero.  Equation  (85) 
cannot  be  established  for  finite  convection  velocity.  However,  as  we  shall  see 
later,  Dyer  (1958)  has  given  a  plausible  argument  leading  to  necessary  conditions 
for  the  validity  of  Eq.  (85)  when  Uc  >  0. 

The  use  of  delta  functions  in  Eq.  (84)  implies  a  flat  wavenumber  spec¬ 
trum  for  wall  pressure  in  the  neighborhood  of  k  =  0.  In  other  -words,  the  plate 
dimension  L  must  be  large  compared  to  all  wavelengths  of  importance  in  the 
driving  wall  pressure  spectrum.  This  circumstance  cannot  hold  if  Phillips' 
result  that  §(k,  in)  -0  as  k  -0  is  accepted.  One  can  see  that  the  obstacles  in 
the  way  of  orderly  analytical  progress  are  indeed  formidable. 

Neglecting  the  effect  of  radiation  loading  on  plate  vibration,  Kraichnan 
found  the  frequercy  spectral  density  of  radiated  sound  power  per  unit  area  to  be 


P(ou)  =  p  c 
o  o 
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Wavenumbers  k>k0  do  not  contribute  to  the  radiated  sound  energy,  for 
they  woi  Id  be  associated  with  a  purely  reactive  radiation  loading  on  the  plate. 
This  fact  is  also  evident  from^Eq.  (63).  One  should  note  that  P(w)  is  not  equiva¬ 
lent  to  the  integral  of  ?f(x2  .  k,  id)  over  k  space.  For  large  xs,  however,  we 
have  the  relation 

k  <kQ  j  ,  2 

7  1  -  (f)  *f(x2,k,  uu)  dk  .  (87) 
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It  would  appear  that  Eq.  (86)  signifies  a  marked  increase  in  the  con¬ 
tribution  to  P(cu)  of  wavenumbers  k  near  the  critical  wavenumber  k0.  However, 
Maidanik  has  shown  (in  an  unpublished  investigation)^  that  the  neglected  radiation 
resistance  becomes  infinite  at  k  =  Icq.  and  hence  V  (k ,  in)  must  become  zero  as 
k  approaches  k0  with  sufficient  rapidity  that  the  contribution  to  P(w)  at  k  =  k0  is 
in  fact  null,  quite  independently  of  the  nature  of  the  driving  wall  pressure 
spectrum  #(k,  in). 

Kraichnan  then  utilized  Eqs. _J82),  (83),  (85),  (86),  and  his  assumed 
form  of  the  wall  pressure  spectrum  §(k,  in)  to  obtain  quantitative  estimates 
of  P(w).  The  estimates  treat  the  contribution  from  turbulence -mean  shear  in¬ 
teraction,  denoted  by  the  approximation  of  Eq.  (36),  independently  of  the  contri¬ 
bution  of  the  turbulence -turbulence  interaction  contributions.  The  latter  are 
found  to  be  important  only  at  low  wavenumbers .  The  concept  of  a  hydrodynamic 
critical  frequency  is  introduced,  i.e.,  the  resonant  frequency  of  the  plate  mode 
with  bending -wave  velocity  equal  to  the  convection  velocity  Uc.  For  the  case  of 
turbulence -mean  shear  driving  wall  pressure  and  critical  structural  damping, 
the  average  radiated  sound  power  per  unit  area, 

CO 

IIr  =  j  P(uu)din 

-  CO 

is  shown  to  be  markedly  greater  when  the  hydrodynamic  critical  frequency  is 
greater  than  the  cutoff  frequency  of  the  driving  spectrum  than  when  it  is  less . 


For  low  constant  damping  and  distributed  convection  velocity,  he  inte  - 
grates  the  radiated  frequency  spectrum  for  the  contribution  from  turbulence -mean 
shear  interaction  between  the  hydrodynamic  critical  frequency  and  the  cutoff  fre¬ 
quency  of  the  driving  spectrum  to  obtain  for  x2  - 00 
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Kraichnan 's  results  show  a  very  strong  dependence  of  far -field  mean 
square  sound  pressure  upon  6*.  This  is  a  direct  consequence^of  the  fact  that  he 
has_assumed  thajt  his  wall  pressure  wavenumber  spectrum  $(k)  -0  as  k  -0.  For 
§  (k)  is  even  in  k  and  hence  (k)  =  O  (k2)  near  k  =  0.  But  §  (k)has  dimensions  of 
pa  x  (length)2  and  hence  §  (k)  ~  p2  5*4  k2  near  k  =  0  .  Kraichnan  assumed  that 
§  ( k)  would  have  this  form  up  to  a  cutoff  wavenumber  k  ^  6* _1  and  be  zero  for 
larger  wavenumbers.  Had  he  chosen  $  (k)  of  the  form  p2  6*  4  k2  exp  (-k2  6*2  ), 
for  example,  he  would  have  obtained  $  (k)  ~p2  b*2  for  k  %  6*"1 .  Since  P  (in)  and 
nr  are  dependent  on  wavenumbers  k  «kp  ,  the  resonant  wavenumbers  of  the 
structure,  it  is  clear  that  the  radiated  sound  is  very  sensitive  to  the  form  of  the 
wall  pressure  spectrum  assumed. 
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A  numerical  example  is  helpful  at  this  point.  Ler  us  consider  a 
simply  supported  structure  of  three -foot -square  steel  plate.,  of  1/4”  thickness 
excited  by  a  boundary  layer  with  Uc  =  20  ft/sec  and  c*  =  0.02  ft.  bor  water 
c0  ?«5000  ft/sec  and  for  steel  C/  17.  000  ft/sec .  The  h vd  rodynamic  critical 
frequency  is 

*h  =  2WcJ  “  2  cps  • 

the  acoustic  critical  frequency  is 

c  2 

f  =  0  °  —  ~  40, 000  cps  , 

B  2  tth  C j- 

the  lowest  plate  resonance  frequency,  corresponding  to'm  ::  n  =  1.  is 

f  1 1  «=  77  cps 

and  the  cutoff  frequency  of  the  wall  pressure  spectrum  is 

U 

f  =  - — t-  %  150  cps  . 

c  .o .  2”  o* 

The  spacing  between  resonant  frequencies  of  the  plate  elements  is 
about  20  cps.  According  to  Kraichnan’s  form  of  i(k),  only  the  lowest  frequencies 
would  contribute  to  P(x)  and  IIr  would_be  proportional  to  ~*4  .  On  the  other  hand, 
the  exponentially  decaying  form  of  5  (k)  would  contribute  to  P(x)  at  higher 
resonant  frequencies  as  well,  with  the  contribution  at  fc>0>  proportional  to  l**  . 

As  a  final  remark  on  Kraichnan’s  work,  we  note  that  he  gave  no 
explicit  consideration  to  the  effect  of  the  decay  of  turbulence  upon  the  wail  pres¬ 
sure  spectrum.  From  physical  considerations,  one  would  expect  that  the  cutoff 
frequency  would  be  significantly  increased  if  there  is  a  finite  decay  time. 

The  response  of  structures  to  a  convecting  field  of  turbulent  boundary 
layer  pressures  has  been  treated  by  an  independent  approach  in  a  related  series 
of  papers:  Lyon  (1956a).  (1956b);  Dyer  (1958),  (1959);  and  Maidanik  (1961). 

In  these  papers  Green’s  function  representations  of  solutions  to  finite  and  infinite 
plate  vibrations  are  used  in  conjunction  with  the  space -time  correlation  repre¬ 
sentation  of  the  driving  wall  pressure.  The  approach  is  thus  distinct  from  the 
Fourier  transform -spectral  density  representations  of  Corcos  and  Liepmann  and 
Kraichnan.  Because  of  the  interrelations  that  exist  between  the  two  approaches, 
results  stemming  from  the  same  physical  assumptions  should  be  comparable. 
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The  Green's  function- "correlation  function  approach  is  characterized 
as  follows: 

(a)  The  wall  pressure  correlation  R  (r,  t  )  is  assumed  to  have 
a  finite  correlation  area  and  can  be  approximated  by  a 
product  involving  delta  functions. 

(b)  The  Green’s  functions  are  expressed  in  infinite  series  of 
the  normal  plate  modes .  Inconsequence,  the  assumption 
is  made  that  contributions  from  different  normal  moues 
are  statistically  uncorrelated. 

(c)  The  effect  of  radiation  loading  on  plate  vibration  is  not 
directly  taken  into  account . 

(d)  Vibrational  displacement  cross  correlations  are  obtained 
that  are  not  averages  over  the  vibrating  surface . 

(e)  The  method  does  not  appear  to  be  particularly  suited  to 
the  calculation  of  an  unbounded  radiation  field  resulting 
from  structural  vibration. 

Items  (b)  and  (c)  are  not  new.  Item  (a)  is  new.  It  cannot  be  logically 
justified  if  Phillips’  theorem  is  accepted.  On  the  other  hand,  its  use  permits  the 
establishment  of  criteria  for  the  validity  of  the  assumption  in  item  (b) .  Item  (d) 
is  advantageous;  it  permits  consideration  of  a  single  structural  clement  rather 
than  an  infinite  array  of  elements  as  used  by  Kraichnan.  Moreover,  one  can 
obtain  the  mean  square  vibratory  displacement  at  a  point  on  the  structure  inst  -ad 
of  only  the  average  of  this  quantity  over  the  surface .  Item  (e),  however,  is  the 
penalty  which  must  be  paid  for  the  advantage  obtained  in  item  (d) .  As  we  shall 
see  later.  Dyer  has  shown  that  item  (e)  is  not  necessarily  a  penalty  in  sonar 
self -noise  problems  where  a  bounded  radiation  field  is  involved . 

Lyon  (1956a)  developed  the  form  of.the  correlation  R  (r,  t  )  from  basic 
principles  of  the  theory  of  probability.  He  considered  a  random  source  field  as 
a  summation  of  "eddies"  centered  in  space -time  with  independent  random  ampli¬ 
tudes,  each  with  the  same  probability  density.  The  probability  density  of 
occurrence  of  an  eddy  in  space  was  assumed  to  be  uniform  over  the  space  in 
question,  and  the  temporal  probability  of  occurrence  of  an  eddy  was  assumed  to 
satisfy  a  Poisson  distribution.  For  the  special  case  where  the  amplitude  of  a 
one-dimensional  eddy  takes  the  values  ±  1  with  equal  likelihood  and  the  eddy 
travels  with  velocity  v  in  the  positive  x  direction  while  decaying  exponentially 
in  time  with  time  constant  9 ,  he  found 

R(r’  T)=i¥  6(r-vT)e‘  (89) 
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where  m  is  the  average  number  of  eddies  created  per  unit  time  per  unit  length . 

He  thus  gave  a  fundamental  theoretical  basis  for  the  limiting  form  we  have  fre¬ 
quently  used  for  the  correlation  function  for  turbulent  wall  pressure  and  provided 
a  probabilistic  framework  for  the  construction  of  more  general  analytical  expres¬ 
sions  for  correlation  functions  to  correspond  to  experimental  results. 

He  pointed  out  further  that  the  response  of  any  linear  system  to  such 
a  random  forcing  function  could,  in  principle,  be  obtained  in  the  following  manner: 

Suppose  0(r,  t)  is  the  response  to  f(r,  t),  i.e. 

*0(r,  t)  =  f<r,  t)  (90) 

where  £  is  a  linear  differential  operator.  Then  one  can  write 


1 

=  / dt<) Jd ?o g <?, t: ;e,to)  f(7,to) 


where  G(r,  t;  rQ,t  )  is  the  Green's  function  of  the  system  satisfying 


<£G(r,t;  r  t  )  =  6  (r-r  )  6(t-t  ) 
o  o  o  o 


and  the  boundary  conditions  of  the  problem .  The  correlation  function  of  the 
response  is  therefore 

t  t’ 


®<?,t)0<?,  f>  dto  yv  jdvo  jd?o 


-  CO  -  00 


xc-(r,t;  ro,to)G(r’,  t';  r;,t^)  f(rQ,t)  f(r;,  t;>  .  (93) 

Lyon  (1956b)  applied  this  methodology  to  the  problem  of  the  excitation 
of  a  string  of  length  L  by  a  moving  field  of  turbulence  with  pressure  correlation 

« 

R  (r,  t)  =  p2  X  6  (r-U  T)e  ‘ ! /3  (94) 


where  X  is  the  correlation  length.  The  linear  operator  in  this  case  is «Sf(r|), 
where  r»is  the  normal  displacement  of  the  string.  The  governing  differential 
equation  is 

Z(r\)  =  h  +  Bq  -  c2q  =  <9 

Ll  L  XX  t'  v 
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Lyon  used  the  Green's  function  in  the  form  of  a  series  of  normal  modes 
of  the  string.  Because  of  the  finite  convection  velocity  Uc,  he  had  to  assume  that 
the  normal  modes  were  statistically  uncorrelated  in  determining  the  modal  mean 
square  velocity  response  of  the  string.  The  principal  result  of  this  investigation 
was  the  discovery  of  the  aforementioned  hydrodynamic  coincidence  effect. 


Lyon's  result  isjjerhaps  more  simply  visualized  by  deriving  the  mean 
square  velocity  response  v2,  where  v  =  r|t,  for  the  case  of  an  infinite  string  by 
applying  the  spectral  approach  to  Eqs.  (94)  and  (95).  One  finds  that 


vz 


uj)  dk  d  iu 


p2  X  2  3  /*  ui2  dk  d  uj 

p2  (2tt)3  J  C(k3c2-ui2)2  +  u>282  ]  [1 +(kU  -  a1)2  B2  ] 

v  c 


(96) 


The  first  bracket  in  the  denominator  of  Eo .  (96)  is  familiar .  It  simply 
indicates  that  the  resonant  response  when  k  =  ^  is  limited  only  by  the  damping 
factor  3.  The  second  bracket,  however,  shows  the  hydrodynamic  coincidence 
effect.  It  has  a  minimal  value  when  k  =  w  .  By  rewriting  this  bracket  in  the 
form  Uc 


2 


one  can  see  that  the  coincidence  effect  is  large  if  c  3>>k~1 .  Thus  one  may  expect 
the  coincidence  effect  to  be  marked  only  for  small  wavelengths .  If  one  turns  now 
tc  the  case  of  a  finite  string,  it  is  intuitively  evident,  and  Lyon  demonstrated  it 
to  be  true,  that  the  modal  mean  square  displacement  will  be  large  for  a  high  mode 

number  n  when  k  =  — —  . 

n  U 

c 

Lyon  made  an  experimental  verification  of  this  result  by  exciting  a  steel 
ribbon  under  tension  with  turbulent  flow  and  measuring  the  reflection  of  a  colli¬ 
gated  beam  of  light  from  the  vibrating  ribbon.  The  reflected  light  actuated  a 
nbototube  whose  output  voltage  was  passed  through  a  1/3  octave  filter,  a  squaring 
circuit  end  integrator,  and  then  to  a  recorder.  The  recorded  signal  was  con¬ 
sidered  to  be  proportional  to  the  mean  square  displacement  of  the  ribbon  for  the 
mode  lying  in  the  1/3  octave  band  selected.  A  strong  hydrodynamic  coincidence 
effect  was  found  for  mode  numbers  n  between  5  and  20,  although  considerable 
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variation  in  Uc  about  the  expected  value  was  found,  depending  upon  inode  number. 
Here  Lyon  assumed  Uc  was  the  mean  stream  flow  velocity.  Recently  McCormick 
(1965)  has  utilized  the  same  basic  technique  as  a  measure  of  Uc  and  reported 
that  he  found  the  ratio  of  Uc  to  mean  stream  velocity  to  be  about  0.8.  This  is  in 
general  conformity  with  the  results  of  Harrison  and  Willmarth. 


Dyer  (1958),  (1959)  extended  Lyon’s  methods  to  the  case  of  a  finite 
damped  plate  excited  by  a  turbulent  boundary  layer  with  a  pressure  correlation 
function  initially  chosen  to  be  of  the  form 

R  (?,  t)  =  5=  exp  f-  i  J  (r,  -  Uc  r)2  +  rf  -  ~-j  (97) 

The  quantity  X  is  the  radius  of  a  correlation  area  for  wall  pressure. 
The  quantity  9  has  the  same  meaning  as  before .  It  is  the  correlation  time  for 
the  moving -axis  autocorrelation  function.  The  frequency  spectral  density  of 
wall  pressure  corresponding  to  Eq.  (97)  can  readily  be  shown  to  be 


where 


(98) 


(99) 


Equation  (98)  can  be  made  to  correspond  reasonably  well  with 
Willmarth ’s  experimental  data  by  a  judicious  choice  of  the  parameter  a>0  . 

Dyer  specialized  the  form  of  his  pressure  correlation  to 

R(r,  T)  =  pp  a  6(ri-Ucr)  6(r3)  e‘-J~L  (100) 

in  order  to  render  his  analysis  tractable . 

The  Green's  function  for  a  finite  plate  is  formally  similar  to  that  for 
a  finite  string,  except  for  the  appearance  of  a  double  series  in  the  normal 
modes  qmn  and  natural  frequencies  wmn.  In  Dyer’s  approach  the  modal  expon¬ 
ential  time  decay  factor  contains  the  effects  of  hysteretic  damping  as  well  as 
Rayleigh  damping.  Hysteretic  damping  is  introduced  by  writing  the  flexural 
rigidity  of  the  plate  as  D  (1  e),  where  e  is  the  hysteretic  loss  factor.  Assum¬ 

ing  statistical  independence  of  modes.  Dyer  calculated  the  modal  normal  velocity 
correlation  as  well  as  the  modal  mean  square  response  of  the  plate.  He  endeavored 
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lu  justify  the  assumption  of  statistical  independence  of  modes  by  a  study  of  the 
double  series  representation  for  the  correlation  function  of  plate  response  where 
the  driving  wall  pressure  correlation  was  of  the  form  of  Eq.  (100).  He  obtained 
a  necessary  condition  for  statistical  independence  of  modes 

U  e  <  <  L  .  (101) 

c 

It  is  evident  that  if  one  chooses  a  wall  pressure  correlation  of  the  form 
of  Eq.  (100),  an  additional  necessary  condition  is 

X  <  <  L  .  *  (102) 


Maidanik  (1961)  studied  the  res  *Mise  of  a  finite  string  to  a  convecting 
random  pressure  field  with  a  space -time  v»rrelation  of  the  form  of  Eq.  (97) 
with  r3  =  0 .  He  obtained  an  explicit  expression  for  the  normal  displacement 
cross -correlation  as  a  doubly  infinite  sum  of  modal  contributions  and  then  showed 
that  the  coupling  terms  are  negligible  compared  with  the  individual  mode  terms  if 
inequalities  (101)  and  (102)  are  satisfied.  He  further  showed  the  delta  function  to 
be  a  good  approximation  to  an  exponential  spatial  correlation  for  the  forcing 
function  for  all  modes  with  wavelengths  much  greater  than  a  correlation  length . 
His  results  obviously  could  be  generalized  to  the  case  of  a  simply  supported  plate 
and  thus  substantiate  Dyer’s  conditions  for  the  statistical  independence  of  modes. 

Dyer’s  results  differ  from  those  of  Kraichnan  in  that  he  emphasized 
the  role  of  the  decay  of  turbulence  by  use  of  a  finite  decay  parameter  9  and  a 
fixed  convection  velocity  Uc .  He  found  again  that  the  modal  mean  square  velocity 
response  of  the  plate  is  maximized  at  hydrodynamic  coincidei  ce.  In  the  case  of 
a  square  plate,  this  occurs  when  the  modal  bending  wave  speed  Cp  given  by 
Eqs.  (64)  and  (78)  is  related  to  the  convection  velocity  by 


(103) 


provided  that  the  inequality  mrr  >  > 


is  satisfied. 


The  inequality  is  physically 


equivalent  to  that  of  Lyon  for  the  case  of  string  excitation  and  requires  that  the 
convected  turbulent  pressure  field  must  decay  in  a  distance  greater  than  an 
m-th  modal  wavelength  if  hydrodynamic  coincidence  is  to  be  important. 


Dyer  pointed  out,  however,  that  in  underwater  applications,  the  fre¬ 
quency  range  of  interest,  the  structural  form  of  sonar  domes,  and  the  convection 
velocities  normally  obtained  conspire  to  make  the  convection  velocity  much  smaller 
than  that  required  for  the  coincidence  effect  to  become  important.  He  therefore 
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directed  his  attention  to  the  case  where 


m  n  U 

— : -  <  <«i>  and  obtained  approxi- 

L  mn 

mate  expressions  for  the  modal  mean  square  response  for  special  cases  of  low 
damping  and  high  frequencies . 


Dyer  continued  his  analysis  to  consider  radiation  into  a  closed  fluid  - 
filled  box.  One  side  of  the  box  was  the  vibrating  plate  and  the  other  sides  were 
considered  to  be  pressure  release  surfaces.  This  case  is  a  reasonable  approxi¬ 
mation  to  the  problem  of  sonar  self  noise  and  is  fortunately  amenable  to  analytic 
solution.  For  the  case  where  the  bending  wave  speeds  of  interest  are  much  less 
than  the  speed  of  sound  in  the  medium,  Dyer  found  that  the  effect  of  interior 
radiation  loading  could  be  included  by  adding  a  modal  mass  impedance  to  that  of 
the  plate.  The  radiation  to  an  outer  un confined  space  involved  the  difficulties  we 
have  noted  erxrlier  and  could  be  treated  only  approximately. 


By  considering  the  average  number  of  modes 


N’ 


V3  L 2  Af 

cxh 


(104) 


in  a  frequency  band  Af  where  is  the  longitudinal  bulk  velocity  in  the  plate. 
Dyer  found  that  the  mean  square  interior  acoustic  pressure  in  a  band  of  width  Af 
was 


(P2) 


band 


U3  V2  Af  . 
mn  mn 


(105) 


Here  A  is  the  correlation  area  and  mn  denotes  a  mode  within  the  band .  Equation 
(105)  refers  to  the  interior  field  pressure  at  the  plate  surface;  hence,  for 
Af  —  0,  P2  /  Af  corresponds  to  5^(0;  k ,  ui)  integrated  over  k  space  rather  than  to 
the  quantity  P(w)  calculated  by  Kraichnan. 

The  work  of  Dyer  and  Maidamk  is  open  to  the  criticism  that  the  use  of 
delta  functions,  or  functions  that  have  the  delta  function  as  a  limiting  form,  to 
represent  the  spatial  correlation  of  wall  pressure  contradicts  Phillips’  theorem. 
We  note,  however,  that  the  only  experimental  data  supporting  Phillips'  theorem 
are  those  of  Bull  .which  show  negative  longitudinal  correlation  R  (rj ,  0,  0)  for 
ra  >  4  6*  .  But  one  could  use  the  model 


„  ,  „  .  —  sin  kt, 

R  (ra,  0,  0)  =  p2  -g— — ^ 

as  an  approximation  to  Bull’s  measurements.  This  function  exhibits  negative 
correlation  and  also  has  the  delta  function  as  its  limiting  form  as  K  -  °  .  It 
appears  that  Maidanik’s  analysis  can  be  adapted  to  this  function  without  altering 
his  conclusions. 
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The  question  of  whether  or  not  Phillips’  theorem  is  satisfied  by  the 
particular  correlation  model  used  in  the  analysis  of  structural  response  does 
not  appear  to  be  especially  important,  if  we  are  interested  in  determining  the 
structural  response  to  a  particular  pressure  field  whose  statistical  properties 
have  been  determined  experimentally.  The  difficulty  comes  when  one  wishes 
to  relate  structural  response  to  properties  that  characterize  a  wide  class  of 
boundary  layer  flows.  Then  the  predicted  dependence  of  structural  response  on, 
say,  the  displacement  boundary  layer  thickness  6*  is  strongly  influenced  by  the 
choice  of  correlation  model . 

Dyer’s  work  is  a  significant  attack  on  the  problem  of  sonar  self  noise, 
but  its  limitations  must  be  recognized .  The  noise  received  by  a  transducer 
located  at  a  distance  from  the  plate  can  be  treated  only  in  an  approximate  manner . 
Furthermore,  the  effect  of  transducer  size  and  directivity  on  the  noise  it  receives 
remains  to  be  studied . 


Ffowcs  Williams  and  Lyon  (1963)  used  the  spectral  approach  to  study 
the  response  of  infinite  and  periodically  simply  supported  plates  to  turbulent 
wall  pressures.  Utilizing  Kraichnan’s  expression  for  the  radiated  power  per  unit 
area 


P(ou) 


V  (k,  uu)dk 

H)T 

L  '  O'  J 


(86) 


they  expressed  V(k,  in)  for  the  infinite  plate  in  terms  of  the  wavenumber -frequency 
spectral  density  of  the  wall  pressure  *(k ,  in)  and  the  plate  admittance  as 


V  (k,  in) 


_ in2  &  (k,  in) _ 

(Pph  ^cpjk4-^  (1  -f-i  £)  j 2 


where  £  is  the  plate  loss  factor . 

$  (k,  in)  = 


They  assumed 

*i(ki)  SbOcs)  *M(m-Ucki) 


(106) 


(107) 


and  inves'igated  the  case  of  low  speed  flow  Uc  ~  0.  They  used  the  experimentally 
derived  relations  for  (kx)  and  $3  (k3  )  determined  from  Rx  (ra)  and  R3  (r3)  of 
Figures  -5  and  4,  respectively: 
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*1  (kx)  = 


k? 


n 


(i+kf 


*3  (ks)  =  i; 


71  a+kixi) 


(108) 


with  X*  ^5.4  6*  and  X3  ^2  6*  based  on  Hodgson's  (1962)  measurements  of  longi¬ 
tudinal  and  lateral  spatial  correlations  of  wall  pressure. 


They  noted  that  k„  >  >  kQ  for  most  low -speed  flow  applications, 
example,  for  a  i/4"  steel  plate  in  water,  kp  >  k0  for  frequencies  below 
fa  40  kcps .  Hence  for  k  <  <  k0  , 


For 


k4  -  k4  (1+is) 

P  1 


ks  , 
P 


/ 


k<  k„ 


§(k)  dk  ~  50  6*  p2  k4 
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from  Eqs.  (108),  and,  consequently: 


P(a>)  —  50 


^  c 

o  o 


&*4  -3 

r  p 
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(109) 


The  numerical  factor  in  Eq .  (109)  is  about  l/20th  that  found  by  Ffowcs 
Williams  and  Lyon .  It  appears  that  an  error  in  calculation  crept  into  their  work . 
However,  the  basic  result  that  the  radiated  power  spectrum  is  proportional  to 
times  the  moving -axis  frequency  spectrum  of  wall  pressure  is  valid.  This  is 
an  extension  of  Ribner's  result,  in  that  far-field  radiation  from  an  infinite  plate 
occurs  only  through  the  mechanism  of  the  decay  of  turbulence.  For,  in  the  case 
of  no  decay,  (w)  =  o(u>)  .  Equation  (109)  then  indicates  only  a  d-c  response. 
This,  of  course,  is  a  consequence  of  setting  Uc  =  0,  for  now  there  is  no  time- 
varying  quantity  left  to  produce  either  a  vibratory  response  or  an  acoustic  radia¬ 
tion. 


For  an  infinite  plate,  Eq.  (109)  indicates  that  P(uj)  is  proportional  to 
6*’,  as  predicted  by  Kraichnan  for  periodically  simply  supported  plates.  This 
follows  from  the  assumption  that  only  the  low- wavenumber  components  of  the  wall 
pressure  contribute  to  far-field  radiation.  These,  in  turn,  are  governed  by 
Eqs.  (108),  which  satisfy  Phillips’  theorem  in  the  sense  that 
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$(k)  -  0  as  k  -*  0  . 


The  situation  for  a  periodically  stiffened  plate  is  markedly  different 
from  that  for  an  infinite  plate.  For  e  <  <  1,  V  (k,  <n)  is  governed  by_resonant 
response  of  the  various  plate  modes;  hence,  we  are  interested  in  $(k,  («)  at 
relatively  high  wavenumbers  r -‘Her  than  at  very  low  wavenumbers.  The  fre¬ 
quency  range  of  interest  for  s-  ’pplications  is 


<  (JU  <  0)  = 

a 


For  this  frequency  range,  hydrodynamic  coincidence  effects  are  not 
important,  and  the  radiated  power  per  unit  area  is  associated  both  with  the  decay 
of  turbulence  and  the  presence  of  plate  boundaries. 


We  can  place  the  analysis  of  Ffowcs  Williams  and  Lyon  in  the  frame¬ 
work  of  that  of  Kraichnan  given  in  Eqs .  (77)  through  (87) .  Neglecting  the  effect 
of  radiation  loading  upon  vibratory  response  and  assuming  the  statistical  inde¬ 
pendence  of  modes,  the  radiated  power  spectral  density  per  unit  area  can  be 
written  as 
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where  the  modal  radiation  efficiency  is 
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The  technique  used  to  evaluate  P(w)  is  similar  to  that  used  by  Dyer. 

One  considers  a  narrowband  Aw  about  wQ  and  obtains 

Aw 

*o+T 

°  P(u>)di»«4 P  C  N  (jrTThjF  Tr)  *1  £  /M4|l»  I*  ^.^RAC 

•'Aw  'op'  >  m.n*' 


w  -  — 
o  2 


inAw 


(113) 


where  the  loss  factor  e  <  <  1  and  the  number  of  resonant  modes  in  Aw  is 


N  = 

4n^ 


(114) 


for  sufficiently  high  frequency  to  .  The  factor 


equal  to 
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w  (p  h)2  2  € 
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Y  (w)|2  dw 
mn  I 


for  any  natural  frequency  of  the  plate  in  the  Aw  band . 

'Die  remainder  of  the  calculation  involves  taking  the  average  of  the 
product  of  two  integrals  over  all  natural  frequencies  in  the  Aw  band .  Both  integrals 
involve  the  dimensionless  functions  jl,  j2  and  jls  j  2  .  Smith  and  Lyon  (1964) 
found 

■  i2  ~  f  m  tt  m  tt  1  m  tt 

h j  »  2l7  [0(kl  "TT)  +  ^  +  T7>J forkl  near±  T7 


(“==)  for  k,j  << 

!  i  mn  s  L,j 


,  m  r  1  or  2  . 


(115) 


!  Is  |  has  a  similar  representation,  with  n  replacing  m  and  L3  replacing  Lx  . 

Equation  (114)  indicates  a  uniform  distribution  of  modes  in  wavenumber 
space  k~  k3,  since  Aw  =  2k  aC,?  Ak.  Hence  the  summation  in  Eq.  (113)  can  be  re¬ 
placed  by  an  integration  over  a  circular  band  in  k,  ks  space.  The  situation 
pertinent  to  the  sonar  self  noise  problem  is  shown  in  Figure  7.  For  wh  <t%  we 
are  dealing  with  modes  that  are  acoustically  slow  but  hyd roaynamica  1 1  y  fast. 

These  modes  include  both  comer  modes  and  edge  modes .  One  can  show,  using 
Eq.  (lie),  tha* 


FIG.  7  CLASSIFICATION  OF  PLATE  MODES  IN 

WAVENUMBER  SPACE 
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corner  modes 


edge  modes 
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where 


In  a  typical  sonar  application,  the  edge -mode  contribution  dominates 
over  that  of  corner  modes.  Ffowcs  Williams  and  Lyon  used  an  average  radiation 
efficiency  CraD  based  upon  edge-mode  contributions  alone.  Maidanik  (1962)  has 
given  more  accurate  calculations  of  modal  radiation  efficiencies. 

ForRM(T)  =  e  'lTl  /0  ,  -  i  .  Equations  (107) 

and  (115)  then  yield 

fN’M2  • 

4  A  " 


An  average  radiation  efficiency  having  been  calculated,  the  summation 
in  Eq.  (113)  can  now  be  replaced  by  an  integration  over  wavenumbers 
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Here  At  (kp)  is  termed  an  effective  correlation  area  of  wall  pressure  by  Ffowcs 
Williams  and  Lyon. 


By  combining  Eqs.  (113)  through  (117),  one  obtains  the  result 
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where 
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K  n  ii  v  C 

'  P*  1 


In  principle,  an  analytic  expression  for  At  (kp)  can  be  obtained  by  per¬ 
forming  the  integration  in  Eq.  (117)  using  the  expressions  for  (kj)  and  $3  (k3) 
given  in  Eq«  (108)  .  Ffowcs  Williams  and  Lyon  chose,  however,  to  perform  a 
numerical  integration  based  upon  Hodgson's  (1962)  data  for  Rt  (r^)  and  the  form 
of  *3  (k3)  given  in  Eq.  (108).  Their  result  is  shown  schematically  in  Figure  8 
together  with  the _At  (kp)  that  would  result  from  Kraichnan’s  (1957)  assumption 
of  the  form  of  $  (k)  . 

If  an  average  value  of  At  proportional  to  is  assumed,  then  in  a 
scaling  situation  where  ojq  9  is  held  constant,  it  follows  from  Eq.  (119)  that 
P(uj0)  ~U3‘5  6*  3/  2  h*3  2  .  The  integrated  sound  power  per  unit  area  will  be 
Jp  (u>0)  d  u)  ~  U  4*5  6*3/2  h*3  2  over  a  fixed  range  of  u>o0  .  Ffowcs  Williams  and 
Lyon  inadvertently  neglected  a  factor  9  in  their  analysis  at  this  point  and  hence 
obtained  proportionalities  increased  by  a  factor  U  6*  from  these. 

An  alternative  interpretation  of  Eq.  (119)  can  be  made  by  applying 
it  and  Figure  8  to  calculate  the  total  acoustic  power  radiated  per  unit  area  for 
a  particular  sonar  application.  For  Uc  =  20  ft/sec,  h  =  1/4",  C^=17,000  ft/sec 
for  steel,  co=5000  ft/sec  for  water  and  6*  =0.2',  we  have  xC^^lOO  sec'1  and 
S  rs0.3  sec.  For  a  frequency  range  of  500  to  2000  cps,  kp6*  lies  between 
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1.12  and  2,24,  From  Figure  -S  wc  may  use  the  approximation 


A  /?>**  **  2  (k  6*)4 
t  P 


If  we  integrate  F.q.  (119)  over  this  frequency  range,  using  p£»;10  °PqU4 
and  assume  L,  =  La  =4  feet,  say.  the  power  radiated  per  unit  area  is  about 


OJ- 

5 


2  /  P  (us)  d  u  «  1.8  x  10’9  watts/cm‘ 


We  note  that  for  Lj  =  L3  =4  feet,  the  frequency  of  the  fundamental 
plate  mode  is  128  rad/sec.  This  is  sufficiently  low  for  the  mode  density  in  the 
frequency  range  u>,  <  tu  <  uu2  to  satisfy  the  conditions  of  the  analysis. 

Assuming  for  simplicity  that  the  flow  parameters  are  uniform  over  a 
sonar  dome  and  that  all  dome  surfaces  are  equidistant  from  an  internal  trans¬ 
ducer,  the  corresponding  sound  pressure  level  of  flow  noise  at  the  transducer 
would  be  36  db  (re  1  pb).  This  prediction  is  much  too  high.  A  more  reasonable 
estimate  can  be  obtained  by  accounting  for  the  effect  of  radiation  loading  upon 
plate  vibratory  levels  in  an  approximate  fashion.  If  in  Eq.  (118)  we  assume  that 
essentially  all  of  the  plate  loss  factor  s  comes  from  radiation  damping,  we  should 
have  £  « 2  ^RAI>  a  quantity  about  50  times  as  great  as  the  value  we  have  used. 

If  this  is  done,  the  sound  pressure  level  at  the  transducer  is  reduced  by  17  db. 
Although  such  an  approach  is  physically  plausible,  it  only  serves  to  emphasize 
that  the  neglect  of  radiation  loading  in  the  analysis  of  plate  vibrations  is  not  a 
good  approximation  when  the  plate  is  water -loaded. 

Two  other  directions  for  improvement  in  the  calculation  of  radiated 
power  per  unit  area  can  be  pointed  out.  The  first  is  analytic .  We  note  that  the 
average  radiation  efficiency  OraD  *s  an  average  over  wavenumbers  involving 
both  edge  and  comer  modes.  Tne  correct  average  can  be  carried  out  directly 
from  Eq .  (1 13)  and  the  total  contributions  from  edge  and  corner  modes  computed . 
The  second  direction  for  improvement  is  essentially  experimental.  From  Fig¬ 
ure  8  and  our  example  above,  it  is  clear  that  the  high -wavenumber  components 
of  §  (k)  are  important  to  the  sonar  problem .  These  components  are  governed  by 
the  spatial  correlations  of  Figures  3  and  4  at  small  values  of  rj/6*  and  r3/6*. 
Experimental  data  at  these  values,  however,  are  lacking.  We  need  essentially 
the  curvatures  of  Rj  (rj)  and  R3  (r3)  at  their  origins. 
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As  a  final  remark,  it  is  now  clear  why  Kraichnan's  predictions  of 
the  radiated  power  per  unit  area  arc  dependent  upon  6*"*  .  Freni  Figure  8  we 
can  see  that  he  had  used  a  form  of  *(k)  proportional  to  6*'*  kp  throughout  its 
entire  non -zero  range.  Such  a  form  would  be  useful  for  calculations  involving 
an  infinite  plate  but  not  a  plate  with  a  practical  resonant  mode  distribution. 
Kraichnan's  analysis,  oi  course,  was  carried  out  before  any  experimental 
data  on  wall  pressure  fluctuations  were  available. 
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Via.  MEASUREMENTS  OF  SOUND  RADIATED  BY  BOUNDARY  - 
LAYER -EX CITED  PANELS 


In  recent  years,  experimental  programs  have  been  conducted  at  sev¬ 
eral  activities  to  measure  both  the  vibratory  response  of  panels  and  the  related 
radiated  sound  power.  Typically  a  thin  panel  is  installed  flush  in  the  wall  of  a 
wind  tunnel  or  air  duct.  A  turbulent  boundary  layer  is  established  on  one  side 
of  the  panel  by  operating  the  tunnel  or  duct  at  subsonic  flow  speeds .  The  vibra  - 
tory  displacement  of  the  panel  is  measured.  In  some  of  the  experiments  the  side 
of  the  panel  away  from  the  flow  is  made  a  portion  of  the  wall  of  a  reverberant 
chamber.  If  the  chamber  is  independently  calibrated  using  a  known  sound  source, 
it  is  possible  to  obtain  the  sound  power  radiated  by  the  panel  by  a  sound  pressure 
measurement  within  the  chamber.  The  sound  power  measured  in  this  manner 
is  the  average  over  the  surface  of  the  panel . 

Characteristically,  these  experimental  programs  encompass  essen¬ 
tially  all  of  the  experimental  measurements  discussed  earlier:  mean  velocity 
profiles,  r.m.s.  wall  pressures,  wall  pressure  frequency  spectra,  wall  pres¬ 
sure  space -time  correlations,  etc.  The  investigators  have  also  endeavored  to 
relate  their  findings  to  one  or  more  of  the  analytical  predictions  considered  here. 

We  shall  discuss  the  results  of  three  series  of  investigations.  The 
first  was  conducted  at  the  University  of  Minnesota  and  reported  by  Tack,  Smith, 
and  Lambert  (1961)  and  Tack  and  Lambert  (1962).  The  second  was  conducted 
at  the  University  of  Toronto  and  reported  by  Ludwig  (1962),  el  Baroudi,  Ludwig, 
and  Ribner  (1963),  and  el  Baroudi  (1964).  The  third  was  conducted  at  the  Air¬ 
plane  Division  of  the  Boeing  Company  and  reported  by  Maestrello  (1965a), 

(1965b),  (1965c),  (1965d).  The  first  investigation  concerned  itself  with  vibra¬ 
tory  response  alone;  the  second  and  third  investigations  were  concerned  with 
both  vibratory  response  and  radiated  sound . 

Since  these  investigations  were  conducted  in  air,  none  are  entirely 
applicable  to  the  sonar  self -noise  problem.  Fluid  loading  is  negligible,  and 
the  hydrodynamic  critical  frequencies  are  significantly  greater  than  the  lowest 
panel  resonant  frequencies,  in  contrast  to  the  usual  situation  in  a  sonar  appli¬ 
cation  .  This  latter  circumstance  is  caused  by  the  great  difficulty  in  obtaining 
a  usable  radiated -sound -to -ambient -noise  ratio  at  very  low  air  flow  speeds  or 
with  thick  panels.  However,  to  the  extent  that  these  experiments  verify  the 
basic  assumptions  and  results  of  the  analytical  studies  of  structural  response, 
they  are  useful  in  substantiating  predictions  for  conditions  pertinent  to  the  sonar 
self  -noise  problem . 
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The  Minnesota  studies  were  based  on  Dyer's  model  of  the  wall  pres¬ 
sure  correlation  function  R(r  ,  s).  Tack,  Smith  and  Lambert  (1961)  measured 
the  properties  of  this  function  and  found  that  Uc  U,  the  centerline  velocity  of 
the  tunnel,  but  that  the  correlation  lengths  >.  i  ,  X3  and  the  decay  time  u.  were 
functions  of  U  and  of  frequency  x  .  Tack  and  Lambert  (1962)  computed  the  root 

mean  square  modal  displacements  J  in  a  manner  equivalent  to  that  of  Dyer 

but  used  a  spectral  rather  than  a  Green's  function  approach.  Working  with  a 

simply  supported  bar,  they  compared  their  calculated  values  of  ./ri^  with 

measured  values  for  mode  numbers  m  =  5,  7,  9,  and  11  and  tunnel  speeds  U 
from  20  to  100  meters /sec.  No  significant  variation  in  displacement  boundary 
layer  thickness  was  possible  in  their  test  set-up. 


The  forms  of  the  plots  of  Jv?  versus  U  for  computations  and  measure- 

m 

ments  were  quite  similar.  However,  it  was  necessary  to  use  v?p  2  =  10  2  q 
to  get  a  good  fit  to  the  experimental  data;  the  more  generally  accepted  value  of 

Jp  2  x  10'°  q  from  Willmarth's  data  gave  calculated  values  of  ,/r(2  that 

were  65%  too  low.  Tack  and  Lambert  credited  this  discrepancy  to  uncertainties 

in  measurements  of  ./p2  .  This  is  certainly  poss-;  e.  A  rough  wall  alone  would 
be  enough  to  account  for  the  difference .  It  is  also  possible  that  the  assumed 
form  of  R(r,  ~)  is  at  fault,  or  that  one  cannot  really  consider  modal  coupling 

to  be  negligible.  We  note,  however,  one  other  possibility:  the  measurements 

of  ./q2~were  made  at  the  center  of  the  bar,  an  anti -modal  point  for  the  modes 

in  question .  The  calculated  value  of  J  T]^  is  an  average  over  the  length  of  the 
bar  and  hence  should  be  significantly  less  than  the  measured  value. 

No  strong  peaking  of  modal  response  was  found  at  hydrodynamic  coin¬ 
cidence  speed,  in  contrast  to  Lyon's  earlier  experimental  results  .  Tack  and 
Lambert  conclude  that  this  lack  of  peaking  is  caused  by  X  \  ,  X  3  and  6  being 
dependent  on  U  and  ju  .  We  recall  that  Kraichnan  achieved  the  same  effect  by 
requiring  Uc  to  be  variable. 

The  experiments  at  the  University  of  Toronto  were  conducted  using 
11" -square  steel  plates  of  thicknesses  .0015",  .002",  .004"  and  .008".  Each 
plate  was  clamped'  in  a  duct  test  section  wall  surrounded  by  a  reverberant 
chamber.  Two  duct  sizes,  12"  by  8"  and  12"  by  1",  were  used  to  permit  the 
effect  of  variation  in  6*  to  be  studied.  Centerline  flow  velocities  up  to  180 ft/sec 
were  used.  Ludwig  (1962)  reported  principally  on  mean  velocity,  wall  pres  - 
sure  and  sound  field  measurements .  El  Baroudi  (1964)  reported  on  wall  pres¬ 
sure  and  plate  vibratory  displacement  space -time  correlation  measurements. 
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The  wall  pressure  spectral  density  measurements  for  the  8"  duct  were 
taken  with  a  transducer  diameter  to  displacement  boundary  layer  thickness  rath. 
d/5*  =  .092.  Dimensionless  plots  of  the  spectral  density  show  a  gradual  down¬ 
ward  slope  up  to  f  6*/U  =  2.5  with  a  sharp  drop  thereafter.  For  the  1"  duct, 
d/6*  =  0.74;  however,  no  clearly  defined  cutoff  frequency  was  found,  although 
the  downward  slope  did  increase  above  f  6*/U  =  0.3.  The  minimal  transducer 
separation  for  both  longitudinal  and  lateral  space -time  wall  pressure  correla¬ 
tion  measurements  was  about  2  6*  ;  hence,  these  data  are  of  limited  utility. 

Radiated  sound  power  measurements  in  the  reverberation  chamber 
showed  that 


P(iu)d, 


h* 


U' 


where  5=5  for  the  .0015"  and  .002"  panels,  5  =  5.5  for  the  .004"  panel  and 
6  <  S  <  8  for  the  .008"  panel.  P(a)  exhibited  an  upper  cutoff  frequency  pro¬ 
portional  to  U/h  for  the  thinner  panels,  but  proportional  to  U2/h  for  the  .008" 
panel. 


Ffowcs  Williams  and  Lyon  (1963)  interpreted  these  experimental 
results  in  terms  of  their  analytical  studies .  They  considered  that  the  principal 
contribution  should  have  come  from  modes  with  frequencies  below  the  hydro - 
dynamical  critical  frequency,  which  for  the  .002”  panel  at  U  =  100  ft /sec  came 
at  1000  cps .  They  concluded  that 


j  P (ju)  d'jj~ 


c*2U5 

h.s 


in  this  circumstance. 


A  calculation  based  on  the  assumptions  that  the  principal  contribution 
came  from  frequencies  above  the  hydrodynamical  critical  frequency  and  that 
A  ~  6*  2  would  yield 


n 

i 

O 


P(uu)  d  ui~ 


U°  6* 

h3/3 


a  result  more  closely  in  accord  with  Ludwig's  experiments.  If,  in  fact,  the 
frequency  range  considered  by  Ludwig  in  his  power  calculations  varied  from 
test  to  test  --as  seems  implied  by  the  signal -to -noise  difficulties  he  encountered, 
dependence  on  U  and  6*  to  fractional  powers  could  result.  Further  clarification 
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of  this  matter  must  await  the  results  of  more  experiments  of  this  type.  The 
weak  dependence  of  radiated  sound  power  upon  6*  found  by  Ludwig  is  somewhat 
questionable.  It  is  based  essentially  on  only  two  values  of  l* .  It  would  have 
been  desirable  for  Ludwig  to  have  related  his  mean  velocity  profiles  to  Coles’ 
law  of  the  wake,  especially  for  the  1"  duct.  This  narrow -duct  flow  may  not 
truly  correspond  to  a  turbulent  boundary  layer  flow . 

El  Baroudi  measured  longitudinal  space -time  and  lateral  spatial 
vibratory  displacement  correlations  and  compared  the  longitudinal  results  with 
correlations  calculated  by  Dyer's  theory.  Isocorrelation  curves  plotted  as 
functions  of  rj  and  Uc~  showed  general  qualitative  agreement  between  measure¬ 
ment  and  calculation .  A  rough  tendency  for  ridges  to  appear  along  the  line 
ri  =  Uc?  could  be  inferred  in  both  cases  for  small  values  of  rx  and  t.  El 
Baroudi  conjectured  that  the  restriction  Uc9  «L  invoked  by  Dyer  to  justify  the 
neglect  of  mode  coupling  could  be  given  the  physical  interpretation  that  vibratory 
waves  running  at  convection  velocity  existed  in  the  plate.  Since  these  forced 
waves  would  decay  at  the  same  rate  as  the  pressure  correlation,  the  restriction 
served  to  insure  that  reflections  would  be  important  only  near  panel  edges .  As 
noted  above,  the  measured  correlations  tend  to  confirm  this  observation  in  a 
qualitative  sense. 

El  Baroudi's  calculations  of  rros  modal  displacement  response  for  a 
simply  supported  plate  included  the  effect  of  modal  coupling.  He  compared  the 
measured  rms  displacement  with  the  total  rms  displacement  calculated  for  all 
modes  considered  to  contribute  significantly.  Comparisons  were  made  at  the 
plate  center  and  at  a  position  near  the  upstream  edge  of  the  plate .  Two  general 
results  were  obtained:  1)  the  contribution  of  modal  coupling  to  the  calculation 
was  small,  and  2)  the  calculated  values  were  consistently  higher  than  the 
measured  values .  The  latter  result  was  attributed  to  Dyer's  use  of  a  delta - 
function  form  for  the  spatial  wall  pressure  correlation,  a  point  source  loading 
being  considered  to  yield  a  greater  response  than  a  distributed  loading. 

Maestrello  carried  out  a  series  of  experiments  at  Boeing  similar  to 
those  conducted  at  the  University  of  Toronto .  He  used  thicker  panels  (0 .020", 

0 .040" ,  0 .060"  and  0 .080")  and  flow  speeds  up  to  700  ft/sec .  His  panels  were 
7"  by  12”  areas  milled  out  of  3/4"  aluminum  plate,  and  die  edge  constraint  was 
more  nearly  cantilever  than  in  the  Tc  ronto  tests .  His  test  setup  did  not  permit 
a  study  of  variation  in  6  * . 

Maestrello  (1965a)  (1965c)  found  statistical  characteristics  of  the 
wall  pressure  essentially  the  same  as  those  of  Willmarth.  He  fitted  them  to  a 
correlation  model  for  R(r,  t)  similar  to  that  of  Dyer,  Eq.  (97).  He  did  not, 
however,  carry  cut  an  analytical  prediction  of  the  vibratory  displacement  or  of 
the  radiated  sound  field. 
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For  free -st ream  Mach  numbers  M  <0.4  he  found  by  reverberation 
chamber  measurements  that 

jp(u))dui~Msh"1,p. 

For  M  >0.4,  however,  he  found 

j  P(yj)d*~M2-3  h'1 

except  for  the  thickest  panel,  which  continued  to  show  M"*  dependence.  He 
attributed  the  change  in  Mach  number  dependence  to  the  fact  that  the  thinner 
panels  exhibited  significant  hydrodynamic  coincidence  effects  at  the  higher  flow 
speeds . 


Maestrello  found  that  P(  uu)  had  both  a  lower  and  an  upper  cutoff  fre  - 
quency.  The  radiated  power  spectrum  fell  off  sharply  outside  of  the  range  set 
by  these  two  frequencies .  The  lower  cutoff  frequency  was  set  by  the  funda¬ 
mental  panel  mode.  The  upper  cutoff  frequency  was  roughly  proportional  to 
U1* 2  h”1*3  .  These  results  appear,  in  general,  to  extend  those  of  Ludwig 
to  panels  of  greater  thickness . 

Maestrello  also  measured  the  longitudinal  space -time  correlation  of 
vibratory  panel  displacement.  His  results  agree  with  those  of  El  Baroudi,  in 
that  isocorrelation  lines  tended  to  be  directed  parallel  to  rx  =  UcT,  indicating 
the  presence  of  forced  vibrational  waves  moving  at  the  convection  speed . 
Correlation  plots  for  the  thickest  panel  tested  did  not  show  this  tendency. 
Maestrello  considered  this  as  further  confirmation  that  the  hydrodynamic 
coincidence  effect  was  not  significant  in  this  case .  His  longitudinal  vibration 
correlation  was  an  even  function  of  t  for  rx  >  0,  implying  the  possibility  of 
forced  vibrational  waves  moving  at  -LL ,  presumably  resulting  from  reflections 
at  panel  edges .  This  interpretation  must  be  considered  as  tentative,  for  one 
could  propose  the  alternative  argument  that  edge  reflections  result  in  resonant 
standing  waves .  Maestrello  is  continuing  his  investigation  of  this  point  by 
measuring  coiTelations  for  panels  of  various  lengths . 

Bull,  Wilby,  and  Blackman  (1963)  determined  experimentally  the  fre¬ 
quency  spectra  of  the  normal  velocity  of  plate  vibration.  They  were  able  to  vary 
c*  independently  of  stream  velocity  U  by  placing  the  test  plates  at  various  dis¬ 
tances  along  their  wind  tunnel  wall.  A  range  0.0525"  <  6*  <0.1775"  was 
obtained  in  this  manner. 

As  we  might  expect,  the  vibration  spectra  they  obtained  showed  strong 
peaks  at  the  plate  resonant  frequencies .  They  found  the  very  interesting  resuit 
that  these  peaks  were,  in  general,  lowered  when  the  boundary  layer  thickness 
was  increased.  The  type  of  analysis  we  conducted  in  Section  VII  would  have  shown 
that,  at  least  for  the  hydrodynamicaily  fast  modes,  the  opposite  should  have 
resulted . 
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IX.  THE  EFFECT  OF  STRUCTURAL  DAMPING  UPON 
PLATE  VIBRATION  AND  RADIATED  SOUND 


We  have  deferred  consideration  of  the  effect  of  structural  damping  on 
plate  vibration  and  the  radiated  sound  field  until  this  point  in  order  to  give  a 
fairly  self-contained  treatment  of  the  matter.  We  saw  earlier  that  structural 
damping  can  be  included  in  the  plate  equation  in  the  form  of  hysteretic  damping 
by  use  of  a  complex  Young’s  modulus  E(l-ie),  with  s  a  material  loss  factor. 

The  presence  of  the  imaginary  i  indicates  that  either  a  discrete  frequency 
excitation  or  a  excitation  frequency  component  of  the  form  e-^1  is  being  con¬ 
sidered.  The  use  of  a  loss  factor  implies  that  stress  in  the  plate  is  a  linear 
function  of  both  strain  and  strain  rate.  The  loss  factor  s  is  found  experimentally 
to  be  a  rather  complex  function  of  the  excitation  frequency »u.  For  this  reason 
modal  loss  factors  are  often  introduced  in  analytical  studies  of  plate  vibration. 
From  analogy  with  the  case  of  a  simple  oscillator,  a  number  of  measures  of 
damping  equivalent  to  the  loss  factor  are  found  in  the  literature,  e.g.,  the 

C  c  s 

percent  of  critical  damping,  100  x  —  where  —  =  ^ ,  and  the  quality,  or  Q, 

1  c  c 

factor  where  Q  =  —  . 

With  less  physical  justification,  a  normal-velocity-dependent,  or 
Rayleigh,  damping  can  be  introduced  by  adding  a  term  of  the  form  Bqt  to  the 
plate  equation.  Although  this  form  of  damping  has  ~>een  used  to  simulate 
structural  damping,  it  is  perhaps  better  suited  to  the  task  of  approximating 
plate  damping  resulting  from  radiation  of  sound  generated  by  the  vibration  itself. 
We  have  already  noted  that  no  general  solution  for  the  problem  of  a  vibrating 
plate  coupled  to  an  unbounded  acoustic  field  has  been  found;  hence,  all  present 
treatments  of  radiation  damping  are  at  best  only  approximations. 

Structural  damping  has  been  included  in  all  of  the  analytical  studies 
of  plate  response  we  have  reviewed  for  the  very  practical  reason  that  if  was 
impossible  to  get  any  answers  without  including  it.  Since  all  of  the  investigators 
essentially  neglected  the  effect  of  radiation  loading,  all  resonant  responses  v'ould 
have  been  infinite  in  the  absence  of  structural  damping. 

Modem  welded  steel  structures  have  very  little  inherent  structural 
damping.  They  tend  to  have  Q’s  of  the  order  of  100  or  more  over  a  wide  fre¬ 
quency  range.  It  is  possible  to  reduce  the  Q's  to  values  of  the  order  of  10,  how¬ 
ever,,  by  use  of  applied  damping  treatments.  The  problem  here  is  to  determine 
what  effect  this  increased  structural  damping  will  have  on  p’  ite  vibration,  and, 
more  pertinently,  on  the  resulting  radiation. 
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Dyer  (1958)  (1959)  was  the  first  to  establish  an  analytical  criterion 
for  effectiveness  of  structural  damping  in  reducing  boundary  layer  noise.  From 
his  study  of  modal  mean  square  vibratory  displacements  lie  found  that  for  fre¬ 
quencies  substantially  greater  than  the  hydrodynamic  critical  frequency  (but  less 
tlian  the  acoustic  critical  frequency)  a  transition  frequency  exists  below  which  a 
given  increase  in  modal  loss  factor  will  cause  a  significant  decrease  in  modal 


mean  square  displacement  q 


mn 


Above  this  frequency  the  same  increase  in  £ 


will  have  little  effect.  Tnis  transition  frequency  is  given  by  Dyer  as 


f 

t 


1 

ncQ 


where  3  is  the  temporal  decay  factor  of  the  turbulent  boundary  layer  and  is 
roughly 


This  transition  frequency  is  not  very  precisely  established.  Tack  and  Lambert 
(1962)  estimated  it  to  be  about  one-half  of  Dyer’s  value.  It  serves,  however,  to 
give  a  quantitative  idea  of  the  frequency  range  in  which  applied  damping  should 
be  an  effective  vibration  reduction  measure.  Dyer  chose  an  example  typifying 
a  sonar  application  where  for  U  =  20  ft/sec  and  6*  =  0.02  ft,  he  found 
9  «  3  x  10"  sec  and  ft  1000  cps.  He  considered  that  an  applied  damping 
treatment  should  be  effective  for  frequencies  from  that  of  the  lowest  plate  mode 
up  to  ft  «  1000  cps  in  this  case. 

For  the  particular  interior  radiation  problem  treated  by  Dyer  (a  rec¬ 
tangular  box  with  all  sides  pressure  release  surfaces  except  the  plate  side), 
the  modal  mean  square  field  pressure  at  the  plate  is  proportional  to  q  ^  . 

Hence  the  transition  frequency  should  be  a  criterion  for  sound  reduction  as  well 
as  for  vibration  reduction.  One  must  be  cautious  about  generalizing  this  near¬ 
field  result;  if  there  is  far -field  radiation  then  there  is  radiation  damping.  It 
may  be  difficult  to  structurally  damp  a  mode  that  already  has  high  radiation 
damping.  On  the  other  hand,  it  is  unnecessary  to  damp  a  mode  with  low  radia¬ 
tion  damping.  Structural  damping  is  most  effective  for  modes  that  lie  between 
these  extremes. 

For  very  low  damping,  either  structural  or  radiation,  the  vibratory 
displacement  or  field  pressure  frequency  spectra  above  the  hydrodynamic 
critical  frequency  are  formed  essentially  from  resonant-mode  responses.  How¬ 
ever,  if  the  damping  is  high,  non -re sonant -mode  responses  become  significant. 
Kraichnan  (1957)  treated  analytically  the  rather  extreme  case  of  critical 


76 


Arthur  23.3Uttlc.Jfnr. 


i 


structural  damping  under  the  assumption  of  statistical  independence  of  modes. 

In  principle,  if  the  damping  is  sufficiently  high  the  modal  resonant  frequencies 
are  shifted  sufficiently  far  from  their  undamped  natural  frequencies  so  that 
modal  coupling  will  result  from  this  reason  alone.  Very  limited  experimental 
evidence  exists  regarding  the  value  of  structural  damping  in  reducing  boundary- 
layer-excited  vibration  and  sound,  but  what  does  exist  is  quite  favorable.  Tack 

and  Lambert  measured  on  a  damped  and  undamped  simply  supported  bar 

over  a  range  of  flow  speeds  from  20  to  100  meters/sec.  For  modes  m  =  5  and 

7  a  reduction  in  Q  from  over  100  to  less  than  10  produced  reductions  in  q  ® 

by  factors  of  from  5  to  10.  The  corresponding  mode  frequencies  fj  «800  cps, 
i-j  «  1500  cps  were  below  ft  «  3000  cps  because  of  the  relatively  high  speeds  and 
thin  boundary  layers  used.  Hence  Dyer’s  criterion  was  not  tested.  If  we 
(cautiously)  translate  these  results  to  a  sound  pressure  level  reduction  for  Dyer’s 
internal  radiation  problem,  we  have  reductions  of  from  14  to  20  decibels  for 
these  modes. 

Maestrello  (1965a)  placed  two  layers  of  damping  tape  on  clamped 
plates  of  0.020”  and  0.080”  thickness  and  measured  the  effect  on  the  total 
radiated  sound  power  level  and  power  spectral  density  in  a  reverberant  cham¬ 
ber.  Damping  applied  to  the  0.020"  panel  reduced  the  sound  power  level  by 
about  18  db,  with  the  principal  reductions  coming  from  frequencies  above 
1000  cps.  The  0.080”  panel  showed  a  sound  power  level  improvement  of  only 
about  3  db,  largely  because  the  response  at  the  fundamental  panel  mode  domin¬ 
ated  the  spectra  and  was  unaffected  by  the  damping  treatment.  Responses  above 
this  fundamental  were  reduced  by  roughly  6  db. 

Dyer’s  transition-frequency  criterion  could  not  be  tested  directly 
against  Maestrello's  results,  as  the  modal  loss  factors  of  the  panels  have  not 
yet  been  reported.  Using  Maestrello's  reported  value  of  9  «  10"*3  sec  and 
estimating  e  «  0.2  we  obtain  ft  «  1500  cps.  However,  very  significant  power 
level  reductions  were  obtained  above  this  value  of  ft- 

It  is  clear  that  more  investigations  of  the  type  conducted  by  Tack  and 
Lambert  and  Maestrello  are  needed  to  clarify  the  role  of  structural  damping  in 
reducing  boundary  layer  noise.  It  would  be  particularly  desirable  from  the  view¬ 
point  of  the  sonar  self  noise  problem  to  conduct  such  tests  in  water  rather  than 
in  air,  so  that  a  more  realistic  assessment  of  the  effect  of  radiation  damping 
could  be  obtained. 
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X.  SUMMARY  AND  CONCLUSIONS 


In  spite  ol'  a  number  of  questions  that  remain  regarding  both  analytical 
and  experimental  results,  substantial  progress  has  b?en  made  in  the  past  dozen 
years  in  understanding  the  basic  mechanisms  of  boundary  layer  noise.  A  num¬ 
ber  of  qualitative  conclusions  are  now  available  for  application  to  the  design  of 
sonar  domes  and  transducers  to  reduce  the  flow  noise  contribution  to  sonar  self 
noise .  It  is  not  too  much  to  expect  that  we  shall  shortly  be  able  to  reduce  these 
conclusions  to  quantitative  form .  This  will  permit  design  decisions  to  be  made 
on  the  basis  of  over-all  signal -to -noise  ratio  for  desired  operating  speeds. 

It  is  difficult  to  single  out  any  one  piece  of  research  as  the  most  sig¬ 
nificant  contribution  to  our  understanding  of  flow  noise.  As  we  have  pointed  out, 
oversimplifications  and,  on  occasion,  errors  have  occurred  in  both  analytical 
and  experimental  studies .  Investigators  disagree  on  a  number  of  significant 
points .  Such  is  to  be  expected .  The  problem  is  one  of  the  most  difficult  in  the 
whole  area  of  applied  mechanics.  However,  each  investigation  we  have  reviewed 
has  contributed  to  the  general  understanding.  Perhaps  the  most  significant  ad¬ 
vances  were  obtained  at  two  points  of  the  development:  the  first  came  with  the 
experimental  determination  of  wall  pressure  statistics,  the  second  with  the  an¬ 
alyses  of  the  response  of  plates  to  wall  pressure  excitation.  The  former  per¬ 
mitted  investigators  to  move  past  the  difficult  theoretical  problem  of  relating 
wall  pressure  to  velocity  fluctuations .  The  latter  served  to  emphasize  those 
characteristics  of  the  boundary  layer  that  are  important  in  determining  the 
vibratory  response  of  plates  and  the  resultant  sound  field . 

A  limited  amount  of  work  remains  to  be  done  on  the  statistics  of  wall 
pressures  for  well  developed  turbulent  boundary  layers  adjacent  to  smooth,  flat 
walls.  For  example,  the  integral  correlation  lengths  contained  in  the  spatial 
correlation  R(r,  o)  are  not  well  defined.  The  most  recent  plots  of  measured 
isocorrelation  curves  for  this  function,  such  as  those  of  Bull,  Wilby.  and  Black¬ 
man  (1963)  and  Maestrello  (1965a),  are  in  sharp  disagreement.  It  is  particu¬ 
larly  important  that  the  curvature  of  R(T,  o)  be  established  in  the  neighborhood 
of  r  =  0.  More  experimental  determinations  of  the  temporal  decay  term  9 
and  the  mean  square  wall  pressure  would  be  desirable  because  of  their 
importance  in  determining  plate  response .  The  question  of  the  resolution  of 
transducers  of  finite  size  remains  in  an  uncertain  status .  It  must  be  settled  - 
hopefully  through  a  novel  experiment  -  as  it  has  an  important  influence  on  both 
the  determination  of  plate  response  and  the  performance  of  flush  -mounted  trans¬ 
ducers  . 


We  have  seen  that  much  of  the  analytical  and  all  of  the  experimental 
work  involving  plate  response  and  re -radiation  has  been  done  for  air  loading 
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rather  than  for  water  loading.  The  analytical  approach  has  been  to  compute  the 
plate  vibratory  response  to  .vail  pressure  fluctuations,  ignoring  ihe  effect  of 
radiation  loading,  and  then  to  determine  separately  the  radiated  power .  Yet 
even  the  crudest  calculation  of  the  resulting  radiation  into  water  shows  that  its 
influence  on  the  vibratory  level  must  be  very  significant.  One  approach  toward 
removing  this  discrepancy  is  indicated  in  the  development  for  an  infinite  plate 
in  Section  VII  leading  to  Eq.  (63).  The  complexity  evident  there  is  troublesome 
but  not  impossible  to  deal  with . 


It  will  probably  be  much  more  difficult  to  analyze  the  forced  vibration 
of  a  plate  with  a  resonant-mode  distribution  under  the  influence  of  radiation 
loading.  Here  the  criteria  for  statistical  independence  of  modes  must  be  care¬ 
fully  re-examined.  If  radiation  damping  is  sufficiently  large,  it  too  will  induce 
significant  mode  coupling.  Criteria  for  mode  independence  will  therefore  have 
to  include  damping  as  weil  as  the  spatial  and  temporal  scales  and  the  convection 
velocity  of  the  turbulent  wall  pressure .  Digital  computation  techniques  will 
probably  have  to  be  applied  more  extensively  than  heretofore  so  that  the  orders 
of  magnitude  of  coupling  terms  can  be  established. 

The  utility  of  structural  damping  as  a  means  of  reducing  flow  noise 
cannot  be  firmly  established  on  an  analytical  basis  until  the  questions  involving 
radiation  damping  have  been  settled.  Here  again,  criteria  for  the  frequency 
range  over  which  structural  damping  is  effective  (such  as  those  of  Dyer)  must 
be  re-examined.  Additional  criteria  involving  the  required  magnitudes  of  struc¬ 
tural  loss  factors  will  have  to  be  developed . 


The  most  practical  experimental  means  of  studying  flow  noise  in 
water  is  by  use  of  buoyant  or  gravity -pr ope.  led  submerged  bodies .  Acoustic 
wavelengths  in  water  are  too  great,  boundaries  are  too  near,  and  ambient  noise 
levels  are  too  high  for  most  water  tunnels  and  tovring  basins  to  be  of  much  serv¬ 
ice.  Most  buoyant  or  gravity -propelled  bodies  used  in  the  past  were  too  small. 
Uncertainties  related  to  Reynolds '-number  scaling  of  their  boundary  layers  and 
to  structural  scaling  were  unavoidably  introduced.  A  large  -  perhaps  quarter 
scale  -  buoyant  or  gravity -propelled  test  vehicle  appears  to  be  a  much  needed 
experimental  tool  for  studying  flow  noise  on  sonar  domes . 

The  paucity  of  experimental  data  relating  surface  roughness  to  flow 
noise  has  been  alluded  to  frequently  in  this  review.  We  recall  that  the  displace¬ 
ment  boundary  layer  thickness  6*  has  been  shown  to  be  important  -  if  contro¬ 
versial  -  parameter  in  the  determination  of  the  radiated  acoustic  power .  Yet 
Coles'  formulation  of  the  laws  of  the  wall  and  of  the  wake  (discussed  in  Section II) 
indicates  that  6*  is  dependent  upon  the  friction  velocity  v#  or,  equivalently,  the 
wall  shear  stress  tq  .  But  these  latter  quantities,  in  turn,  depend  upon  surface 
roughness.  Clearly,  an  analytical  and  experimental  evaluation  of  the  relation  of 
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surface  roughness  to  radiated  acoustic  power  is  needed  if  one  is  to  answer  the 
practical  question  of  how  smooth  a  sonar  dome  must  be  to  keep  now  noise  within 
acceptable  limits . 


The  effect  of  surface  roughness,  however,  is  only  one  aspect  of  the 
boundary  layer  that  requires  further  study.  Profitable  directions  for  further 
work  on  wall  pressure  statistics  lie  in  the  application  of  experimental  techniques 
already  developed  for  the  study  of  more  general  flows.  A  limited  effort  in  this 
direction  has  been  carried  out  at  Southampton  by  Bull,  Wilby,  and  Blackman (1963), 
who  measured  wall  pressures  in  the  transition  region  between  laminar  and  turbu¬ 
lent  boundary  layers.  Such  studies  should  be  extended  to  deal  with  adverse  pres¬ 
sure  gradients,  flows  in  interference  regions,  and  separated  flows.  Although 
Bull  found  no  significant  change  in  mean  square  pressure  in  the  neighborhood  of 
normal  transition,  the  question  is  open  as  to  whether  or  not  sharp  increases  in 
fluctuating  pressures  will  occur  where  transition  has  been  delayed  by  use  of 
compliant  coatings  or  other  means . 

A  substantial  clarification  of  the  problem  of  structural  response  and 
re-radiation  of  sound  was  achieved  when  the  analytical  model  of  a  periodically 
simply  supported  plate  was  introduced.  We  have  seen  that  for  frequencies  of 
interest  in  sonar  applications  phenomena  related  to  hydrodynamic  or  acoustic 
coincidence  are  quite  unimportant .  The  radiated  acoustic  power  is  produced 
primarily  by  resonant  response  of  structural  modes  to  wall  pressure  fluctuations . 
It  is  sensitive  to  the  distribution  of  resonant  frequencies,  to  radiation  and  struc¬ 
tural  damping,  and  to  certain  statistics  of  the  wall  pressure,  namely  the  wave- 
number  spectrum  at  high  wavenumbers  and  of  the  moving-axis  frequency  spec¬ 
trum  of  the  wail  pressure.  The  customarily  measured  fixed-axis  frequency 
spectrum  is  of  itself  relatively  unimportant,  except  in  relation  to  flush-mounted 
sonar  transducers . 

As  useful  as  the  model  of  a  periodically  simply  supported  plate  has 
been,  we  must  recognize  its  limitations .  Actual  stiffened  sonar  domes  have 
much  more  complex  boundary  conditions .  They  may  have  modes  of  vibration 
which  contribute  significant  radiation  that  cannot  be  modeled  in  this  manner. 

In  mathematical  term0,  eigenfunctions  which  are  products  of  sinusoids  do  not 
form  a  complete  set ,  Further  investigations  involving  more  general  boundary 
conditions  would  be  desirable . 

Finally  there  remains  a  significant  gap  between  the  analyses  and  ex¬ 
periments  yielding  acoustic  radiation  from  plates  and  the  determination  of  how 
this  radiation  affects  the  self  noise  as  measured  by  a  sonar  array.  Two  differ¬ 
ent  aspects  of  this  radiation  have  been  predicted  and  measured,  namely,  the 
sound  pressure  level  and  the  sound  power  level  per  unit  area,  together  with 
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their  corresponding  frequency  spectra,  it  would  appear  that  the  sound  pres¬ 
sure  information  is  useful  primarily  in  predicting  near-field  effects  such  as 
would  occur  when  a  sonar  transducer  is  placed  very  near  the  dome  surface. 
On  the  other  hand,  the  sound  power  information  should  yield  essentially  the 
reverberant  field  level  existing  throughout  a  central  region  of  the  dome .  The 
dome,  however,  is  not  really  a  reverberant  chamber,  and  account  should  also 
be  taken  of  shadowing  and  directivity  effects  associated  with  large  sonar  ar¬ 
rays  that  occupy  a  significant  fraction  of  the  dome  volume .  Analysis  of  these 
aspects  of  the  flow  noise  problem  should  be  directed  toward  determining  the 
significant  measurements  to  be  taken  in  full-scale  experiments  with  sonar 
installations  aboard  ships .  Very  little  guidance  of  this  sort  is  evident  from 
the  current  literature  of  the  problem . 
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u(y)>u\(y2) 
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SECTION  II 


free  stream  velocity  in  positive  x  or  xx  direction 

mean  flow  velocity  in  direction  of  U  parallel  to  a 
wall;  y  =  y2  is  a  distance  normal  to  the  wall 

universal  constants  of  the  law  of  the  wall 


friction  velocity  = 


wall  shear  stress 


V? 
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* 

6 
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density  of  quiescent  fluid 

boundary  layer  thickness 
kinematic  viscosity  of  the  fluid 
thickness  of  the  laminar  sublayer 

displacement  boundary  layer  thickness 
momentum  boundary  layer  thickness 


R 


x 


friction  coefficient  = 


i  s 
!pou 


distance  along  a  wall  in  direction  of  U 


Reynolds  number  = 


Ux 
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*  Listed  in  order  of  appearance  in  the  text. 
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P(u1,u2) 
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^(k) 
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p(x,t) 


R(r,T) 


§{k) 


UJ 

V 
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Coles'  wake  function 

form  nnrnmeter  of  Coles*  lav/  cf  the  wake 

SECTION  m 

displacement  vectors,  x  =  x.,  r  = 

til 

fluctuating  velocity  component  in  i  direction 

spatial  velocity  correlation  tensor  =  u.(x)u.(x+r) 

joint  probability  density  function 
wavenumber  vector 

energy  spectrum  tensor  or  velocity  wavenumber 
spectral  density  tensor 

wavenumber  transform  of  u.(  x  ) 
time 

fluctuating  pressure 
time  variable 

space -time  pressure  correlation  =  p(x,t)p(x  + 
wavenumber  spectral  density  of  pressure 
radial  frequency 

frequency  spectral  density  of  pressure 
frequency  cross  spectral  density  of  pressure 
integration  volume 
integration  time 
convection  velocity 
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R( r, t  ) 
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ccpTtrwT  iw 


uuu  i  ivn  4  v 


fluid  density 

speed  of  sound  in  quiescent  fluid 
Kronecker  delta  (=  1  for  i  =  j,  -  0  for  i  ^  j) 
LighthiU’s  turbulence  tensor 
fluid  stress  tensor 


viscous  components  of  fluid  stress 
surface  area 

retardation  operator  t— *-i  -  r/cQ 

distance,  either  If]  or  i  x  -  y! 
unit  inward  normal  vector  to  surface  S 


stress  vector  exerted  on  the  fluid  by  a  surface  element 

dS  P.  =  n.T.. 
i  j  iJ 

wall  area 


SECTION  V 

fluctuating  velocity  component  normal  to  a  wall 
moving  frame  coordinates,  >T  =  x,  -  U^t 

moving  axis  wall  pressure  cross  correlation  function 
moving  axis  pressure  temporal  autocorrelation  function 
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S(xs ,  xj ,  k,  x’) 


2 

V  P 
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^-f> 

u(rx,f) 
V(rx, 0 

p  (x,t) 
o 

f 

o 

Af 

ait) 


B(t) 


distance  from  boundary  layer  transition  point  in  direc- 

flAM  A  f  f  T 

V  *V/  1  «  VA  V-/ 

2 

d  u.u. 

wavenumber -frequency  spectrai  density  of  P  - — 

o  dx.dXj 


SECTION  VI 


root -mean -square  wall  pressure 
dynamic  pressure  =  |  Pq  Vs 

experimental  frequency  spectral  density  of  wall  pres 

—  A  r»  JVmA 

UXV.  . . 


r.. 


U) 


frequency  in  cps  =  — 


experimental  cross  spectral  density  of  wall  pressure 

real  part  of  §(r x ,  f) 

imaginary  part  of  ?(rx ,  f) 

filtered  output  wall  pressure  signal 

filter  center  frequency 

filter  bandwidth 

component  of  downstream  wall  pressure  contributing 
toUCr^f) 

component  of  downstream  wall  pressure  contributing 
to  V(rx ,  f) 

uncorrelated  component  of  downstream  wall  pressure 
experimental  frequency  spectral  density  of  Y(t) 
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transducer  diameter 


d 

Ro(r,  •  )  fiitered  output  cross  correlation  of  wall  pressure 


a  phase  angle  of  cross  spectral  density  of  wall  piessure 

=  r.t-’VU 
1  c 

C(r,£)  amplitude  function  of  filtered  output  cross  correlation 

of  wall  pressure 

/rix\  _  Ai®\ 

Al-y-l,  Bl-y-  J  Corcos' similarity  functions,  C(r,u»)%  A  I  ~jj~  1  B 
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Hr) 

y(k) 

Ji 

?(k,T) 


r.(x,t) 


D 


transducer  radius 

mean  square  wall  pressure  measured  by  a  transducer 
of  finite  area 

frequency  spectral  density  of  wall  pressure  measured 
by  a  transducer  of  finite  area 

transducer  geometry  function 

wavenumber  transform  of  9(  r  ) 

Bessel  function  of  first  kind  of  order  one 

wavenumber  transform  of  R(r,T  ) 

wavenumber  transform  of  R(  r,  t  )  as  measured  by 
transducers  of  finite  area 

SECTION  vn 

normal  displacement  of  a  plate 

radiated  acoustic  pressure  above  and  below  a  plate 

1  Rh^ 

plate  flexural  rigidity  =  -rr  - 

12  (l-o2) 
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k 

P 


Yniintr#Q  m/vltiluc 

- - O . . ** 

plate  thickness 
Poisson's  ratio 
plate  density 

material  loss  factor 
quiescent  fluid  density 

fluid  velocity  potential 

wavenumber -frequency  spectral  density  of  field  pres¬ 
sure  P 

u) 

acoustic  wavenumber  =  — 

c 
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=  V  k2  -  k2 
o 

real  part  of  Y 

w 

plate  wavenumber  =  — 

c 

P 


c 
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plate  free  bending  wave  speed  =  V 

P 

displacement  =  rl  -  U  t 
correlation  radius  of  wall  pressure 


longitudinal  (bulk)  speed  of  sound  in  plate 
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v(X,X’) 


V  (X) 
mn 


F  (x) 
mn 


Y  (a) 

mn 


V  (*) 
mn 


V(k,x) 


0(ka) 


radius  of  hound  region  of  plate  used  by  Corcos  and 
Liepmann 

space -time  correlation  of  plate  acceleration 
length  of  side  of  a  square  plate 
orthonormal ized  plate  eigenfunction 

plate  eigenfrequency 

frequency  transform  of  plate  velocity 
frequency  transform  of  plate  displacement 
modal  velocity  coefficient 

modal  wall  pressure  coefficient 

modal  plate  admittance 

Rayleigh  damping  coefficient  for  plate 

modal  frequency  spectral  density  of  wall  pressure 

modal  frequency  spectral  density  of  normal  plate 
velocity 

correlation  area  of  wall  pressure 

wavenumber -frequency  spectral  density  of  normal 
plate  velocity 

frequency  spectral  density  of  radiated  sound  power  per 
unit  area 

radiated  sound  power  per  unit  area  averaged  over  the 
plate  and  over  time 

order  of  magnitude  of  k* 
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hydrodynamic  critical  frequency  =  75—; — 

ZTlAC-j 

i  2 

2  ll 

plate  radius  of  gyration  *  =  — 

average  number  of  e  ic  ies  per  unit  time  per  unit  length 
eddy  velocity 

time  constant  for  ;dd y  decay 
linear  differential  operator 
response  function 
forcing  function 
Green's  function  related  to  £. 

wave  velocity  of  a  string  c  =  T/P > 

v 

string  tension 

lineal  density  of  a  string 

natural  frequency  of  a  string 

mode  wavenumber  of  a  string 

average  number  of  plate  modes  ir»  a  frequency  band  Af 

longitudinal  wavenumber  spectral  density  of  wall 
pressure 

longitudinal  correlation  length  of  wall  pressure 
lateral  wavenumber  spectral  density  of  wall  pressure 
lateral  correlation  length  of  wall  pressure 
radial  hydrodynamic  critical  frequency  =  2uf^ 
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radial  acoustic  critical  frequency  =  2  -f 

< 

finite  plate  weighting  functions 


plate  dimensions 

average  number  of  plate  resonant  modes  in  ^ ,  =  2^’ 
modal  radiation  efficiency 

acoustic  wavelength  at  the  acoustic  critical  frequency 

effective  correlation  area  of  wall  pressure  at  plate 

wavenumber  k 
P 

radiation  loss  factor 


driving  point  impedance  of  an  undamped  infinite  plate 

=  3P  hxC; 
p 


damping  parameter  =  £•-*> 
mean  free  bending  wave  path 


2n 


plate  area 
plate  perimeter 


bar  eigenfunction 


Mach  number  =  — 
c 

o 


97 


3rilmr  3nr. 


SECTION  IX 


x  100 


percent  of  critical  damping 


quality  factor  =  - 


Dyer’s  transition  frequency  = 


3itliur 


